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The 1937 theoretical discovery of Majorana fermions — whose defining property is that they are their own 
anti-particles — has since impacted diverse problems ranging from neutrino physics and dark matter searches to 
the fractional quantum Hall effect and superconductivity. Despite this long history the unambiguous observation 
of Majorana fermions nevertheless remains an outstanding goal. This review article highlights recent advances 
in the condensed matter search for Majorana that have led many in the field to believe that this quest may soon 
bear fruit. We begin by introducing in some detail exotic 'topological' one- and two-dimensional superconduc- 
tors that support Majorana fermions at their boundaries and at vortices. We then turn to one of the key insights 
that arose during the past few years; namely, that it is possible to 'engineer' such exotic superconductors in the 
laboratory by forming appropriate heterostructures with ordinary s-wave superconductors. Numerous proposals 
of this type are discussed, based on diverse materials such as topological insulators, conventional semiconduc- 
tors, ferromagnetic metals, and many others. The all-important question of how one experimentally detects 
Majorana fermions in these setups is then addressed. We focus on three classes of measurements that provide 
smoking-gun Majorana signatures: tunneling, Josephson effects, and interferometry. Finally, we discuss the 
most remarkable properties of condensed matter Majorana fermions — the non-Abelian exchange statistics that 
they generate and their associated potential for quantum computation. 



I. INTRODUCTION 

Three quarters of a century ago Ettore Majorana introduced 
into theoretical physics what are now known as 'Majorana 
fermions': particles that, unlike electrons and positrons, con- 
stitute their own antiparticles.^ The monumental significance 
of this development required many intervening decades to 
fully appreciate, and despite being an 'old' idea Majorana 
fermions remain central to diverse problems across modern 
physics. In the high-energy context, Ettore 's original sugges- 
tion that neutrinos may in fact be Majorana fermions endures 
as a serious proposition even today.^ Super symmetric theo- 
ries further postulate that bosonic particles such as photons 
have a corresponding Majorana 'superpartner' that may pro- 
vide one of the keys to the dark matter puzzle.^ Experiments at 
the large hadron collider are well-positioned to critically test 
these predictions in the near future. Condensed matter physi- 
cists, too, are fervently chasing Majorana's vision in a wide 
variety of solid state systems, motivated both by the pursuit 
of exotic fundamental physics and quantum computing ap- 
plications. While a definitive sighting of Majorana fermions 
has yet to be reported in any setting, there is palpable opti- 
mism in the condensed matter community that this may soon 
change. ISI^ 

Unlike the Majorana fermions sought by high-energy 
physicists, those pursued in solid state systems are not fun- 
damental particles — the constituents of condensed matter are, 
inescapably, ordinary electrons and ions. This fact severely 
constrains the likely avenues of success in this search. In con- 
ventional metals, for example, electron and hole excitations 
can annihilate, but since they carry opposite charge are cer- 
tainly not Majorana fermions. In operator language this is 
reflected by the fact that if cj. adds an electron with spin a, 
then its Hermitian conjugate Ccr is a physically distinct oper- 
ator that creates a hole. If Majorana is to surface in the solid 
state it must therefore be in the form of nontrivial emergent 
excitations. 



Superconductors (and other systems where fermions pair 
and condense) provide a natural hunting ground for such ex- 
citations. Indeed, because Cooper pair condensation sponta- 
neously violates charge conservation, quasiparticles in a su- 
perconductor involve superpositions of electrons and holes. 
Unfortunately, however, this is not a sufficient condition for 
the appearance of Majorana fermions. With only exceedingly 
rare exceptions superconductivity arises from 5-wave-paired 
electrons carrying opposite spins; quasiparticle operators then 
(schematically) take the form d = uc\. + vc^, which is still 

physically distinct from = ^*c| + u*c^. Thus whereas 
charge prevents Majorana from emerging in a metal, spin is 
the culprit in conventional 5- wave superconductors. 

As the preceding discussion suggests, 'spinless' 
superconductors — i.e., paired systems with only one active 
fermionic species rather than two — provide ideal platforms 
for Majorana fermions. By Pauli exclusion. Cooper pairing 
in a 'spinless' metal must occur with odd parity, resulting in 
p-wave superconductivity in one dimension (ID) and, in the 
most relevant case for our purposes, p -\- ip superconductivity 
in two dimensions (2D). These superconductors are quite 
special: as Sec. [ll| describes in detail, they realize topological 
phases that support exotic excitations at their boundaries 
and at topological defects.^ Most importantly, zero-energy 
modes localize at the ends of a ID topological p-wave 
superconductoi^, and bind to superconducting vortices in 
the 2D p -\- ip casJ^l. These zero-modes are precisely the 
condensed matter realization of Majorana fermions^ that 
are now being vigorously pursued. 

Let 7 denote the operator corresponding to one of these 
modes (the specific realization is unimportant for now). This 
object is its own 'anti-particle' in the sense that 7 = 7^ and 
7^ = 1. We caution, however, that labeling 7 as a particle — 
emergent or otherwise — is a misnomer because unlike an ordi- 
nary electronic state in a metal there is no meaning to 7 being 
occupied or unoccupied. Rather, 7 should more appropriately 
be viewed as Sifractionalized zero-mode comprising 'half of 



2 



a regular fermion. More precisely, a pair of Majorana zero- 
modes, say 7i and 72, must be combined via / = (71 +i72)/2 
to obtain a fermionic state with a well-defined occupation 
number. While this new operator represents a conventional 
fermion in that it satisfies f ^ P and obeys the usual anti- 
commutation relations, / remains nontrivial in two critical re- 
spects. First, 7i and 72 may localize arbitrarily far apart from 
one another; consequently / encodes highly non-local entan- 
glement. Second, one can empty or fill the non-local state de- 
scribed by / with no energy cost, resulting in a ground-state 
degeneracy. These two properties underpin by far the most in- 
teresting consequence of Majorana fermions — the emergence 
of non-Abelian statistics. 

A brief digression is in order to put this remarkable phe- 
nomenon in proper perspective. Exchange statistics character- 
izes the manner in which many-particle wavefunctions trans- 
form under interchange of indistinguishable particles, and is 
one of the cornerstones of quantum theory. There indeed 
exists a rather direct path from particle statistics to the ex- 
istence of metals, superfluids, superconductors, and many 
other quant uni p hases, not to mention the periodic table as 
we know it.^^^Sl it has long been appreciated that for topo- 
logical reasons 2D systems allow for particles whose statis- 
tics is neither fermionic nor bosonic.^"^ Such anyons come 
in two flavors: Abelian and non-Abelian. Upon exchang- 
ing Abelia n anyons — which arise in most fractional quantum 
Hall stateJI^IiSiMnl — the wavefunction acquires a statistical 
phase e*^ that is intermediate between —1 and 1. Non-Abelian 
anyons are far more exotic (and elusive); under their exchange 
the wavefunction does not simply acquire a phase factor, but 
rather can change to a fundamentally different quantum state. 
As a result subsequent exchanges do not generally commute, 
hence the term 'non-Abelian'. An important step toward find- 
ing experimental realizations of the second flavor came in 
1991 when Moore and Read introduced a set of 'Pfaffian' trial 
wavefunctions for fractional quantum Hall states that support 
non-Abe lian any ons . ^ ^ — Several theoretical and experimen- 
tal works^lS^ indicate that the observed quantum Hall state 
at filling factor^^ = 5/2 may provide the first realization 
of such a non-Abelian phase. Read and Green^^ later pro- 
vided a key breakthrough that in many ways served as a step- 
ping stone for the new directions reviewed here. In particular, 
these authors established an intimate connection between the 
superficially very different Moore-Read Pfaffian states and a 
topological spinless 2D p + zp superconductor — deducing that 
universal properties of the former such as non-Abelian statis- 
tics must also be shared by the latter (which crucially can arise 
in weakly interacting systems). 

With this backdrop let us now describe how non-Abelian 
statistics arises in a 2D spinless p -\- ip superconductor. Con- 
sider a setup with 2N vortices binding Majorana zero-modes 
71,..., 2 AT. One can (arbitrarily) combine pairs of Majoranas 
to define N fermion operators fj = (72^-1 + i72j)/2 cor- 
responding to zero-energy states that can be either filled or 
empty. Thus the vortices generate 2^ degenerate ground 
states^ ^ that can be labeled in terms of occupation numbers 



= fjfj by 

|ni,n2, . . . ,nAr). (1) 

Suppose that one prepares the system into an arbitrary ground 
state and then adiabatically exchanges a pair of vortices. 
Because this process swaps the positions of two Majorana 
modes, each being 'half of a fermion, the system generally 
ends up in a different ground state from which it began. More 
formally the exchange unitarily rotates the wavefunction in- 
side of the ground- state manifold in a non-commutative fash- 
ion. The vortices — because of the Majorana zero- mode s that 
they bind — therefore exhibit non-Abelian statistics. I^*^^^ 

One might naively conclude that in this regard the Majo- 
rana zero-modes bound to the ends of a ID topological p- 
wave superconductor are substantially less interesting than 
those arising in 2D. After all, exchange statistics of any type 
is ill-defined in ID because particles inevitably 'collide' dur- 
ing the course of an exchange.^ This is the root, for instance, 
of the equivalence between hard-core bosons and fermions in 
ID. Fortunately this obstacle can be very simply surmounted 
by fabricating networks of ID superconductors; envision, say, 
an array of wires forming junctions, with topological p-wave 
superconductors binding Majorana zero-modes interspersed 
at various locations. Such networks allow the positions of 
Majorana zero-modes to be meaningfully exchanged, which 
remarkably still gives rise to non-Abelian statistics despite 
the absence of vortices. f^^Jtlllxhus ID and 2D topological su- 
perconductors can both be appropriately described as non- 
Abelian phases of matter. [As an interesting aside, Teo and 
Kane first showed that non-Abelian statistics can even appear 
in three dimensions, where exchange has long been assumed 
to be trivial.^^^ 

The observation of Majorana fermions in condensed mat- 
ter would certainly constitute a landmark achievement from a 
fundamental physics standpoint, both because it could mean 
the first realization of Ettore Majorana's theoretical discovery 
and, far more importantly, because of the non-Abelian statis- 
tics that they harbor. Moreover, success in this search might 
ultimately prove essential to overcoming one of the grand 
challenge s in t he field — the synthesis of a scalable quantum 
computer.E^l^^lESI xhe basic idea is that the occupation num- 
bers rij = 0, 1 specifying the degenerate ground states of 
Eq. Q can be used to encode 'topological qubits'.^^ Cru- 
cially, this quantum information is stored highly non-locally 
due to the arbitrary spatial separation between pairs of Majo- 
rana modes corresponding to a given rij . Suppose now that 
temperature is low compared to the bulk gap; if manipula- 
tions are carried out adiabatically the system then essentially 
remains confined to the ground-state manifold. The user can 
controllably manipulate the state of the qubit by adiabatically 
exchanging the positions of Majorana modes, owing to the ex- 
istence of non-Abelian statistics. In principle the environment 
can also induce (unwanted) exchanges, thereby corrupting the 
qubit, but this happens with extraordinarily low probability 
due to the non-locality of such processes. This is the basis of 
fault-tolerant topological quantum computation schemes that 
elegantly beat decoherence at the hardware level. While 
braiding of Majorana fermions alone permits somewhat lim- 
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ited topological quantum information processing,'^ the addi- 
tional unprotected operations needed for universal qua ntum 
computation come with unusually high error 

thresholds.EMI 

The search for Majorana fermions is thus fueled also by the 
potential for revolutionary technological applications down 
the road. 

In the beginning of this introduction we noted that re- 
searchers are optimistic that this search may soon come to 
fruition. One might reasonably wonder why, given that we 
live in three dimensions, electrons carry spin, and p-wave 
pairing is scarce in nature. To a large extent this opti- 
mism stems from the recent revelation that one can engineer 
low-dimensional topological superconductors by judiciously 
forming heterostructures with conventional bulk 5-wave su- 
perconductors. This new line of attack could eventually lead 
to 'designer topological phases' persisting up to relatively 
high temperatures, perhaps measuring in the lOK range or 
beyond. The conceptual breakthrough here originated with 
the seminal work of Fu and Kane in the context of topologi- 
cal insulators ,E§EOi which paved the way for many subsequent 
proposals of a similar spirit. We devote a large fraction of this 
review — Sees. [Ill] and |IV] — to discussing these new routes to 
Majorana fermions. 'Classic' settings such as the = 5/2 
fractional quantum Hall state and Sr2Ru04 (which of course 
remain highly relevant to the field) will also be discussed, but 
only briefly. An omission that we regret is a discussion of 
Helium-3, where seminal work related to this subject was car- 
ried out early on by Volovik and others; see the excellent book 
in Ref. 8 SectionjV] explores the key question of how one ex- 
perimentally identifies Majorana modes once a suitable topo- 
logical phase is fabricated. The long-term objectives of ob- 
serving non-Abelian statistics and realizing quantum compu- 



conductor. This model's many virtues include the fact that in 
this setting Majorana zero-modes appear in an extremely sim- 
ple and intuitive fashion. Following Kitaev, we introduce op- 
erators Cx describing spinless fermions that hop on an A/'-site 
chain and exhibit long-range-ordered p-wave superconductiv- 
ity. The minimal Hamiltonian describing this setup reads 



H 



-E 

2 ^ 



[tcl.Cx 



+1 



(2) 

where /i is the chemical potential, t > is the nearest- 
neighbor hopping strength, A > is the p-wave pairing am- 
plitude and (j) is the corresponding superconducting phase. For 
simplicity we set the lattice constant to unity. 

It is instructive to first understand the chain's bulk prop- 
erties, which can be conveniently studied by imposing peri- 
odic boundary conditions on the system (thereby wrapping 
the chain into a loop and removing its ends). Upon passing to 
momentum space and introducing a two-component operator 
Cl = [c^, c_/e], one can write H in the standard Bogoliubov- 
de Gennes form: 



H 



E 

keBZ 



\ ^k 



-^k 



(3) 



with e/e = —tcosk — /i the kinetic energy and Aj^ = 
—iAe'^^ sin k the Fourier-transformed pairing potential. The 
Hamiltonian becomes simply 



k^k 



keBZ 



tation are taken up in Sec.lVll Finally, we offer some closing when expressed in terms of quasiparticle operators 



thoughts in Sec. |VII[ For additional perspectives on this fas- 
cinating problem we would like to refer the reader to several 
other reviews and popular articles: Refs. [31471 fT2l [T3l [5TH531 



II. TOY MODELS FOR TOPOLOGICAL 
SUPERCONDUCTORS SUPPORTING MAJORANA MODES 

This section introduces toy models for topological ID and 
2D superconductors that support Majorana fermions. We will 
explore the anatomy of the phases realized in these exotic su- 
perconductors and elucidate how they give rise to Majorana 
modes in some detail. Later parts of this review rely heavily 
on the material discussed here. Indeed, our perspective is that 
all of the recent experimental proposals highlighted in Sees. 



(4) 



(5) 



^k = T^, , Vk 



Ik 



|A| V2E^k ^ V A 

where the bulk excitation energies are given by 



£^buik(^) = V e 



Uk, (6) 



(7) 



III and IV are, in essence, practical realizations of these toy 
models. The ideas developed here will also prove indispens- 
able when we discuss experimental detection schemes in Sec. 
[viand non-Abelian statistics in Sec.lVll 



A. ID spinless p-wave superconductor 

We begin by reviewing Kitaev 's toy lattice modeP, intro- 
duced nearly a decade ago, for a ID spinless p-wave super- 



Equation ^ demonstrates that the chain admits gapless bulk 
excitations only when the chemical potential is fine-tuned to 
fi = tor —t, where the Fermi level respectively coincides with 
the top and bottom of the conduction band as shown in Fig. 
[TJa). The gap closure at these isolated /i values reflects the p- 
wave nature of the pairing required by Pauli exclusion. More 
precisely, since A ^ is an odd function of k, Cooper pairing at 
k = or k = ±7r is prohibited, thereby leaving the system 
gapless at the Fermi level when jj. = ±t. Note that the phases 
that appear at /i < —t and jj. > t slyq related by a particle- 
hole transformation; thus to streamline our discussion we will 
hereafter neglect the latter chemical potential range. 

The physics of the chain is intuitively rather different in the 
gapped regimes with ji < —t and < t — the former con- 
nects smoothly to the trivial vacuum (upon taking ji —oo) 
where no fermions are present, whereas in the latter a partially 
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FIG. 1. (a) Kinetic energy in Kitaev's model for a ID spinless p- 
wave superconductor. The p-wave pairing opens a bulk gap except at 
the chemical potential values /x = ±t displayed above. For > t 
the system forms a non-topological strong pairing phase, while for 
< t a topological weak pairing phase emerges. The topological 
invariant u distinguishing these states can be visualized by consid- 
ering the trajectory that h(/c) [derived from Eq. ( [lo| )] sweeps on the 
unit sphere as k varies from to tt; (b) and (c) illustrate the two types 
of allowed trajectories. 



filled band acquires a gap due to p-wave pairing. One can 
make this distinction more quantitative following Read and 
Green^^ by examining the form of the ground-state wavefunc- 
tion in each regime. Equation ^ implies that the ground state 
|g.s.) must satisfy a/c|g.s.) = for all /c so that no quasiparti- 
cles are present. Equations ([5]) and ^ allow one to explicitly 
write the ground state as follows, 



ig.s.)oc n [i + ^c.p.{k)ciAm 



^C.p.(^) 



0</c<7r 
Vk_ _ 
Uk 



(8) 



where |0) is a state with no Ck fermions present with mo- 
menta in the interval < < tt. One can loosely inter- 
pret (^c.p.(^) as the wavef unction for a Cooper pair formed 
by fermions with momenta k and —k. An important differ- 
ence between the fi < —t and \fi\ < t regimes is manifested 
in the real-space form (fc.pX^) — J^^^^^^CpX^) large 



Iv^c.p.(^)l 



const, 



li< -t 
\^\<t 



(strong pairing) 
(weak pairing) . 



It follows that ji < —t corresponds to a strong pairing 
regime in which 'molecule-like' Cooper pairs form from two 
fermions bound in real space over a length scale whereas 
in the weak pairing regime < t the Cooper pair size is 
infinite^^. We emphasize that this distinction by itself does 
not guarantee that the weak and strong pairing regimes con- 
stitute distinct phases. Indeed, similar physics occurs in the 
well-studied "BCS-BEC crossover" in 5-wave paired systems 
where no sharp transition arises^^^^. The fact that the weak 
and strong pairing regimes are distinct phases separated by 



a phase transition at which the bulk gap closes is rooted in 
topology. 

There are several ways in which one can express the 'topo- 
logical invariant' (akin to an order parameter in the theory of 
conventional phase transitions) distinguishing the weak and 
strong pairing phaseP. We will follow an approach that 



closely parallels the 2D case we address in Sec. II B Let us 
revisit the Hamiltonian in Eq. ([3]), but now allow for addi- 
tional perturbations that preserve translation symmetry. The 
resulting 2x2 matrix 1-Lk can be expressed in terms of a vector 
of Pauli matrices a = a^x + a^y + a^z as follows. 



Hk = Hk) • a 



(10) 



for some vector h{k). (A term proportional to the identity 
can also be added, but will not matter for our purposes.) Al- 
though we are considering a rather general Hamiltonian here, 
the structure of h{k) is not entirely arbitrary. In particu- 
lar, since the two-component operator Ck in Eq. ^ satisfies 



a^Ck, the vector h{k) must obey the important 



relations 



hx,y{k) = -hx^y{-k), h^{k) = h^{-k). 



(11) 



Thus it suffices to specify h{k) only on the interval < k < 
TT, since h(/c) on the other half of the Brillouin zone follows 
from Eq. ([TT]). 

Suppose now that h(/c) is non-zero throughout the Brillouin 
zone so that the chain is fully gapped. One can then always 
define a unit vector h{k) = h{k)/\h{k)\ that provides a map 
from the Brillouin zone to the unit sphere. The relations of 
Eq. ^U) strongly restrict this map at A: = and tt such that 



h(0) 5oz, h(7r) = s^z, 



(12) 



where sq and s^^ represent the sign of the kinetic energy (mea- 
sured relative to the Fermi level) at A: = and tt, respectively. 
Thus as one sweeps k from to tt, h(/c) begins at one pole of 
the unit sphere and either ends up at the same pole (if 50 = ^7^) 
or the opposite pole (if sq = —s^). These topologically dis- 
tinct trajectories, illustrated schematically in Figs. [TJb) and 
(c), are distinguished by the Z2 topological invariant 



(13) 



which can only change sign when the chain's bulk gap closes 
[resulting in h(A:) being ill-defined somewhere in the Brillouin 
zone].^^ Physically, u = +1 if at a given chemical potential 
there exists an even number of pairs of Fermi points, while 
u = —1 otherwise. From this perspective it is clear that u = 
+1 in the (topologically trivial) strong pairing phase while 

= — 1 in the (topologically nontrivial) weak pairing phase. 

The nontrivial topology inherent in the weak pairing phase 
leads to the appearance of Majorana modes in a chain with 
open boundary conditions, which we will now consider. The 
new physics associated with the ends of the chain can be most 
simply accessed by decomposing the spinless fermion opera- 
tors Cx in the original Hamiltonian of Eq. ([2]) in terms of two 
Majorana fermions via 



-i0/2 



-{7b,x -^i7A,x)' 



(14) 
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FIG. 2. Schematic illustration of the Hamiltonian in Eq. |T6| when 
(a) / 0, t = A = and (b) /i = 0, t = A / 0. In the 
former limit Majoranas 'pair up' at the same lattice site, resulting 
in a unique ground state with a gap to all excited states. In the lat- 
ter, Majoranas couple at adjacent lattice sites, leaving two 'unpaired' 
Majorana zero-modes 7^,1 and 7s, at at the ends of the chain. Al- 
though there remains a bulk energy gap in this case, these end- states 
give rise to a two-fold ground state degeneracy. 



The operators on the right-hand side obey the canonical Ma- 
jorana fermion relations 



(15) 



In this basis H becomes 



N 



x=l 
N-1 



- T ^ [(A + t)jB,xlA,x+l + (A - t)jA i(]16) 

X = l 

Generally the parameters /i, t, and A induce relatively com- 
plex couplings between these Majorana modes; however, the 
problem becomes trivial in two limiting cases^. 

The first corresponds to /i < but t = A = 0, where the 
chain resides in the topologically trivial phase. Here the sec- 



ond line of Eq. ([16]) vanishes, leaving a coupling only between 
Majorana modes ^a,x and 7^ at the same lattice site as Fig. 
|2ja) schematically illustrates. In this case there is a unique 
ground state corresponding to the vacuum of fermions. 
Clearly the spectrum is gapped since introducing a spinless 
fermion into the chain costs a finite energy Note that this 
is entirely consistent with our treatment of the chain with pe- 
riodic boundary conditions; in the trivial phase the ends of the 
chain have little effect. We emphasize that these conclusions 
hold even away from this fine-tuned limit provided the gap 
persists so that the chain remains in the same trivial phase. 

The second simplifying limit corresponds to /i = and 
t = A ^ 0, where the topological phase appears. Here the 
Hamiltonian is instead given by 



t ^"^ 

H = -i- 7B,a^7A,x+i, 



(17) 



as Fig.|2jb) illustrates. In terms of new ordinary fermion oper- 
ators dx = ^ (7A,a^+i +^7B,a^), the Hamiltonian can be written 



N-1 



X=l 



(18) 



In this form it is apparent that a bulk gap remains here 
too — consistent with our results with periodic boundary 
conditions — since one must pay an energy t to add a dx 
fermion. However, as Fig. |2jb) illustrates the ends of the 
chain now support 'unpaired' zero-energy Majorana modes 
7i = 7^ 1 and 72 = 7b, at that are explicitly absent from 
the Hamiltonian in Eq. ([17]). These can be combined into an 
ordinary — though highly non-local — fermion. 



/ = 2(71+^72), 



(19) 



which couples Majorana fermions only at adjacent lattice sites 



that costs zero energy and therefore produces a two-fold 
ground- state degeneracy. In particular, if |0) is a ground state 
satisfying /|0) = 0, then |1) = /^|0) is necessarily also a 
ground state (with opposite fermion parity). Note the stark 
difference from conventional gapped superconductors, where 
typically there exists a unique ground state with even parity so 
that all electrons can form Cooper pairs. 

The appearance of localized zero-energy Majorana end- 
states and the associated ground-state degeneracy arise be- 
cause the chain forms a topological phase while the vacuum 
bordering the chain is trivial. (It may be helpful to imag- 
ine adding extra sites to the left and right of the chain, with 
ji < — t for those sites so that the strong pairing phase forms 
there.) These phases cannot be smoothly connected, so the 
gap necessarily closes at the chain's boundaries. Because this 
conclusion has a topological origin it is very general and does 
not rely on the particular fine-tuned limit considered above, 
with one caveat. In the more general situation with /i 7^ 
and t ^ A (but still in the topological phase) the Majorana 
zero-modes 71 and 72 are no longer simply given by 7^,1 and 
7b, at; rather, their wavefunctions decay exponentially into the 
bulk of the chain. The overlap of these wavefunctions results 
in a splitting of the degeneracy between |0) and |1) by an en- 
ergy that scales like e~^l^, where L is the length of the chain 
and ^ is the coherence length (which diverges at the transition 
to the trivial phase). Provided L ^ however, this splitting 
can easily be negligible compared to all relevant energy scales 
in the problem; unless specified otherwise we will assume that 
this is the case and simply refer to the Majorana end- states as 
zero-energy modes despite this exponential splitting. 

Finally we comment on the importance of the fermions be- 
ing spinless in Kitaev's toy model. This property ensures that 
a single zero-energy Majorana mode resides at each end of the 
chain in its topological phase. Suppose that instead spinful 
fermions — initially without spin-orbit interactions — formed a 
p-wave superconductor. In this case spin merely doubles the 
degeneracy for every eigenstate of the Hamiltonian, so that 
when < t each end supports two Majorana zero-modes, 
or equivalently one ordinary fermionic zero-mode. Unless 
special symmetries are present these ordinary fermionic states 
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will move away from zero energy upon including perturba- 
tions such as spin-orbit coupling. (Note that even for a spin- 
less chain it is in principle possible for multiple nearby Majo- 
rana modes to coexist at zero energy if certain symmetries are 
present; see Refs. [5911611 for examples. Time-reversal sym- 
metry can also protect pairs of Majorana end- states in 'class 
Dlir ID superconductors with spin.^^^i^ 

This by no means implies that it is impossible to experimen- 
tally realize Kitaev's toy model and the Majorana modes it 
supports with systems of electrons (which always carry spin). 
Rather these considerations only imply that a prerequisite to 
observing isolated Majorana zero-modes is lifting Kramer's 
degeneracy such that the electron's spin degree of freedom be- 
comes effectively 'frozen out'. We will discuss several ways 
of achieving this, as well as the requisite p- wave superconduc- 



tivity, in Sec. Ill 



B. 2D spinless p + ip superconductor 

In two dimensions, the simplest system that realizes a topo- 
logical phase supporting Majorana fermions is a spinless 2D 
electron gas exhibiting p+zp superconductivity. We will study 
the following model for such a system. 



H 



2m 



+ ^ [e'U{d,^idy)^^H.c.] 



(20) 



where ip^ (r) creates a spinless fermion with effective mass m, 
/i is the chemical potential, and A > determines the p-wave 
pairing amplitude while (j) is the corresponding superconduct- 
ing phase. For the moment we take the superconducting or- 
der parameter to be uniform, though we relax this assumption 
later when discussing vortices. To understand the physics of 
Eq. ( [20| ) we will adopt a similar strategy to that of the previ- 
ous section — first identifying signatures of topological order 
encoded in bulk properties of the p -\- ip superconductor, and 
then turning to consequences of the nontrivial topology for the 
boundaries of the system. 

In a system with periodic boundary conditions along x and 
y (i.e., a superconductor on a torus with no edges) translation 
symmetry allows one to readily diagonalize Eq. ( [20| by going 
to momentum space. Defining ^(k)^ = [^/^^(k), ?/^(— k)],one 
obtains 

with e{k) = 1^ - /i and A(k) = iAe^^(fc^ + iky). A 
canonical transformation of the form a(k) = ii(k)?/;(k) + 
v(k)?/;^ (— k) diagonalizes the remaining 2x2 matrix. In terms 
of these quasiparticle operators the Hamiltonian reads 



^=2 



H 



-I 



(27r)2 



£buik(k)aT(k)a(k). 



(22) 



(a) \ 

topological \ 


/ 

2m / 
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— V >■ 


(weak pairing)V 
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(strong pairing) 
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Ttrivial) 




FIG. 3. (a) Kinetic energy for a spinless 2D electron gas exhibiting 
p + ip superconductivity. The pairing opens a bulk gap except when 
/i = 0. This gapless point marks the transition between a weak 
pairing topological phase at /i > and a trivial strong pairing phase 
at /i < 0. These phases are distinguished by the Chem number C 
which specifies how many times the map h(k) [derived from Eq. 
( [26t ] covers the entire unit sphere as one sweeps over all momenta 
k. As |k| increases from zero, in the trivial phase h(k) covers the 
shaded area in (b) but then 'uncovers' the same area, resulting in 
Chern number C = 0, whereas in the topological phase the map 
covers the entire unit sphere once as illustrated in (c) leading to | C| = 
1. 



The coherence factors ii(k) and v(k) take the same form as 
in Eq. ([6]), and the bulk excitation energies are similarly given 
by 



^buik(k) = Ve(^)' + |A(k) 



(23) 



For any /i > the bulk is fully gapped since here the pairing 
field A(k) is non-zero everywhere along the Fermi surface. 
As one depletes the band the bulk gap decreases and eventu- 
ally closes at /i = 0, where the Fermi level resides precisely 
at the bottom of the band as shown in Fig. [3ja). (The gap 
closure here arises because Pauli exclusion prohibits p-wave 
pairing at k = 0.) Further reducing /i reopens the gap, which 
remains finite for any /i < 0. 

As in the ID case the intuitively different /i > and /i < 
gapped regimes can be quantitatively distinguished by exam- 
ining the ground- state wavefunction^^, which can be written 
as 



|g.s.) oc 



V^c.p.(k) = 



n [1 

kx>0,ky 

u{k) 



.^c.p.(k)^(-k)tV(k)t]|0) 



\ A 



(24) 



where |0) is a state with no V^(k) fermions present with non- 
zero momentum. The 'Cooper pair wavefunction' (pc.p.O^) 
again encodes a key difference between the /i > and a < 
regimes. In real space one finds the asymptotic forms^ 



l^c.p.W 



-|r|/C 

v.l-1 



/i < (strong pairing) 
/i > (weak pairing). 



(25) 



demonstrating that /i < corresponds to a 'BEC-like' strong 
pairing regime, whereas with /i > a 'BCS-like' weakly 
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paired condensate forms from Cooper pairs loosely bound in 
space. 

Also as in the ID case, topology underlies the fact that the 
weak and strong pairing regimes constitute distinct phases that 
cannot be smoothly connected without closing the bulk gap. 
To expose the topological invariant that distinguishes these 
phases, consider a 2D superconductor described by a Hamil- 
tonian of the form^^ 



H(k) = h(k) • a 



(26) 



with h(k) a smooth function that is non-zero for all momenta 
so that the bulk is fully gapped. One can then define a unit vec- 
tor h(k) that maps 2D momentum space onto a unit sphere. 
Assuming that h(k) tends to a unique vector as |k| ^ oo (in- 
dependent of the direction of k), the number of times this map 
covers the entire unit sphere defines an integer topological in- 
variant given formally by the Chern number 

C = J^[h-{dkhxdkh)]- (27) 

The integrand above determines the solid angle (which can 
be positive or negative) that h(k) sweeps on the unit sphere 
over an infinitesimal patch of momentum space centered on 
k. Performing the integral over all k yields an integer that 
remains invariant under smooth deformations of h(k). The 
Chern number can change only when the gap closes, making 
h(k) ill-defined at some momentum. 

Consider now the Hamiltonian in Eq. ( [2T] ) for which 
h^{k) = Re[A(k)], hy{k) = Im[A(k)], and h,{k) = e{k). 
Notice that for momenta with fixed |k|, and hy always 
sweep out a circle on the unit sphere at height hz. As |k| 
increases from zero in the /i < strong pairing phase, hz 
begins at the north pole, descends towards the equator, and 
then returns to the north pole as |k| oo. Thus in the (topo- 
logically trivial) strong pairing phase h(k) initially sweeps 
out the shaded region in the northern hemisphere of Fig. [3]^b) 
but then 'unsweeps' the same area, resulting in a vanishing 
Chern number. In contrast, for the (topologically nontrivial) 
/i > weak pairing phase hz transitions from the south pole 
at k = to the north pole when |k| oo; the map h(k) 
therefore covers the entire unit the sphere exactly one time as 
shown schematically in Fig.[3jc), leading to a nontrivial Chern 
number C = —1. [Note that other integer Chern numbers are 
also possible. For instance, dip — ip superconductor carries a 
Chern number C = +1 in the topological phase. An /-wave 
superconductor with A(k) oc {kx -\- iky)^ provides a more 
nontrivial example. In this case for momenta with fixed |k|, 
hx and hy trace out a circle on the unit sphere three times, 
yielding a Chern number C = —3 in the weak pairing phase 
(see, e.g., Ref. 66).] 

We will now explore the physical consequences of the non- 
trivial Chern number uncovered in the topological weak pair- 
ing phase. Consider the geometry of Fig.|4ja), where a topo- 
logical p-\-ip superconductor occupies the annulus and a triv- 
ial phase forms elsewhere. We will model this geometry by 
H in Eq. (20) with a spatially dependent /i(r) that is positive 



inside the annulus and negative outside. Since these regions 
realize topologically distinct phases one generically expects 
edge states at their interface, which we would like to now un- 
derstand following various authors^^*^^^. Focusing on low- 
energy edge modes and assuming that /i(r) is slowly varying, 
one can discard the — V^/(2m) kinetic term in H. A minimal 
Hamiltonian capturing the edge states can then be written in 
polar coordinates (r, 0) as 



edge 



i/j + H.c. 



(28) 



Because of the e*^ factor above, the p -\- ip pairing field cou- 
ples states with orbital angular momentum quantum numbers 
of different magnitude. In what follows it will be convenient 
to gauge this factor away by defining t/j = e~*^/^?/;^ (Note 
that ide ^ ide -\- 1/2 under this change of variables, though 
the constant shift vanishes in the pairing term by Fermi statis- 
tics.) Crucially, the new field ijj^ must exhibit anti-periodic 
boundary conditions upon encircling the annulus. 

In terms of (r) = [?/;'^(r), ?/^'(r)], the edge Hamiltonian 
becomes 



J^edge = d2r*'t(r)^(r)vl/'(r), 

n{v)-. 



-li{r) Ae-'*(-a, 

To find the edge state wavefunctions satisfying H(r)x(r) 
Ex{t^), it is useful to parametrize x(r) as 

-^/2[/(r)+z^(r)]- 



.(29) 



Xn(r) 



Jc/>/2 



[/(r) - zg{r)] 



(30) 



where n is a half -integer angular momentum quantum number 
to ensure the proper anti-periodic boundary conditions. The 
functions / and g obey 

{E^nA/r)f = -i[^{r)-Adr]g 
{E-nA/r)g = i[^{r)^Adr]f. (31) 
For modes well-localized at the inner/outer annulus edges, it 



suffices to replace r ^ R 



in/out 



on the left-hand side of Eqs. 



(31 ). Within this approximation one finds that the energies of 
the outer edge states are 

nA 



out 



while the corresponding wavefunctions follow from / 
and [/i(r) — Adr\g = 0. The latter equations yield 



°"'(r) = e 



-te 



•412 



(32) 
= 

(33) 



which indeed describes modes exponentially localized around 
the outer edge. Similarly, the inner-edge energies and wave- 
functions are given by 
nA 

^in 



Eu 



(34) 



(35) 
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FIG. 4. (a) A topological p -\- ip superconductor on an annulus sup- 
ports chiral Majorana edge modes at its inner and outer boundaries, 
(b) Energy spectrum versus angular momentum n for the inner (red 
circles) and outer (blue circles) edge states in the setup from (a). Here 
n takes on half-integer values because the Majorana modes exhibit 
anti-periodic boundary conditions on the annulus. An ^ flux pierc- 
ing the central trivial region as in (c) introduces a branch cut (wavy 
line) which, when crossed, leads to a sign change for the Majorana 
edge modes. The flux therefore changes the boundary conditions to 
periodic and shifts n to integer values. This leads to the spectrum 
in (d), which includes Majorana zero-modes 71 and 72 localized at 
the inner and outer edges. The two- vortex setup in (e) supports one 
Majorana zero-mode localized around each puncture, while the outer 
boundary remains gapped. 



Figure |4jb) sketches the energies versus angular momentum n 
for the inner (red circles) and outer (blue circles) edge states. 

These edge modes exhibit several remarkable features. 
First, they are chiral — the inner modes propagate clockwise 
while the outer modes propagate counterclockwise, as is clear 
from Fig. |4jb). (For sl p — ip superconductor, the chiralities 
are reversed.) While this is reminiscent of edge states found 
in the integer quantum Hall effect, there is an important dis- 
tinction. The edge states captured above correspond to chiral 
Majorana modes which, roughly, comprise 'half of an inte- 
ger quantum Hall edge state. To be more precise let us expand 
^'(r) in terms of edge-mode operators Fn 



[x;r(r)r-+xr(r)c 



(36) 



Since the upper and lower components of ^'(r) are related 
by Hermitian conjugation, Eqs. ( 33]|36 ) imply that Fn = 
(F^^{°^*)^ This property in turn implies that (i) only edge 
modes with energy ^ > [solid circles in Fig.|4jb)] are phys- 
ically distinct, and {ii) the real- space operators 

pin/out ^ gin^pin/out ^ jpin/out f ^37^ 



are in fact Majorana fermions. 

While these chiral Majorana edge modes become gap- 
less when the topological and trivial regions are thermody- 
namically large, in any finite system there remains a unique 
ground state. This is a direct consequence of the anti-periodic 
boundary conditions on ^'(r) which led to half-integer val- 
ues of n and hence minimum edge-excitation energies of 
A/(2i?in/out)- The physics changes qualitatively when a flux 
quantum $ = |^ threads the central trivial region as shown 
in Fig.Qc). This flux induces a vortex in the superconducting 
pair field so that (say) A Ae~*^ in Eq. (28). The edge 
Hamiltonian in the presence of this vortex can be written in 
terms of our original fermion fields ^^(r) = [V^^(r), V^(r)] 
(which exhibit periodic boundary conditions) as 



d2r^t(r)H(r)^(r) 



(38) 



with H(r) again given by Eq. (29). The Hamiltonians with 
and without a vortex appear identical, so the edge-state en- 



ergies and wavefunctions again take the form of Eqs. (32 



35 ), with one critical difference. Since ^(r) exhibits periodic 
boundary conditions, the angular momentum quantum num- 
ber n now takes on integer values. The edge state spectrum 
sketched in Fig.|4jd) then includes two zero-energy Majorana 
modes 71 and 72, one localized at each interface. These Ma- 
jorana zero-modes are the counterpart of the Majorana end- 
states discussed in Sec. II A and similarly result in a two-fold 
ground state degeneracy for the p-\-ip superconductor. (Tech- 
nically, the edge- state wavefunctions overlap if the topologi- 
cal region is finite, splitting this ground-state degeneracy by an 
energy that is exponentially small in the width of the annulus. 
Throughout we will neglect such a splitting unless specified 
otherwise.) 

The shift in boundary conditions underlying the formation 
of Majorana zero-modes can be intuitively understood as fol- 
lows. First, note that sending ijj e'^^^^'^tp is equivalent to 
changing the phase of the superconducting pair field by Scj). 
Thus a (5^ = 27r shift in the superconducting phase, while ir- 
relevant for Cooper pairs, effectively leads to a sign change for 
unpaired fermions [such as the edge mod e op erators Fn ' 



see the wavefunctions in Eqs. (33) and ([35])]. To account 
for such sign changes it is useful to take the superconducting 
phase in the interval [0, 27r) and introduce branch cuts indicat- 
ing where the phase jumps by 27r. The wavy line in Fig.Qc), 
for instance, represents the branch cut arising due to the |^ 
flux. A Majorana fermion crossing that branch cut acquires a 
minus sign, thereby changing the anti-periodic boundary con- 
ditions to periodic as we found above in our analytic solution. 
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This perspective is exceedingly valuable partly because it al- 
lows one to immediately deduce where Majorana zero-modes 
form even when an analytic treatment is unavailable. In the 
two- vortex setup of Fig. |4je), for example, chiral Majorana 
edge states at the inner boundaries exhibit periodic boundary 
conditions and therefore host zero-modes, whereas the outer 
edge modes suffer anti-periodic boundary conditions and ex- 
hibit a finite-size gap. Furthermore, this picture will prove es- 
sential for understanding interferometry experiments and non- 
Abelian statistics later in this review. 

So far we have discussed chiral Majorana modes resid- 
ing at fixed boundaries of a topological p -\- ip superconduc- 
tor. An interface between topological and trivial regions can 
also form dynamically when a magnetic flux penetrates the 
bulk of a (type II) topological superconductor. In this case 
the vortex core — which has a size of order the coherence 
length ^ ^ VF/ikp^), with vp and hp the Fermi velocity 
and momentum — forms the trivial region. Adapting Eq. ( [34| ) 
to this situation, the energies of the chiral Majorana modes 



bound to an |^ vortex are given roughly by 



\n\{kF/\f 



(39) 



where kp^ is the bulk gap, Ep is the Fermi energy, and n 
takes on integer values.^^ The spectrum of Eq. ( 39 ) reflects the 
p-\-ip analog^ ^ of Caroli-de Gennes-Matricon states^ ^ bound 
to vortices in 5- wave superconductors.^^ Since n is an integer 
the vortex binds a single Majorana zero-mode (unlike the 5- 
wave case where all bound states have finite energy). It is im- 
portant to observe, however, that this zero-mode is separated 
by a 'mini-gap' £^mini-gap ~ {kp^Y / Ep from the next ex- 
cited state. In a 'typical' superconductor £^mini-gap can easily 
be a thousand times smaller than the bulk gap, which can pose 
challenges for some of the proposals we will review later on. 
In this regard, an appealing feature of the ID p-wave super- 



conductor discussed in Sec. II A is that there the Majorana 
zero-modes are generally separated from excited states by an 
energy comparable to the bulk gap. 

Because we considered a spinless p ^ ip superconductor 
above, each |^ vortex threading a topological region binds 
a single localized Majorana zero-mode. Remarkably, stable 
isolated Majorana zero-modes can also form in a spinful 2D 
electron system exhibiting spin-triplet p-\-ip superconductiv- 
ity. For such a superconductor the pairing term in Eq. pQ| 
generalizes to^ 



triplet 



/ 



d'r- 



H.c. 



^^CTy{d-a){d:,^idy)^- 

(40) 

where V^J, (r) creates an electron with spin a =t, | and spin 
indices are implicitly summed. Note that i^tripiet is invari- 
ant under arbitrary spin rotations about the d direction, i/j 

d o-^^ but transforms nontrivially under all other spin ro- 
tations, reflecting the spin-triplet nature of Cooper pairs. In 
the presence of an ordinary |^ vortex, the superconducting 
phase (j) rotates by 27r around the vortex core. This vortex 
binds a pair of Majorana zero-modes (one for each electron 
spin) which generically hybridize and move to finite energy 



upon including spin-mixing perturbations such as spin-orbit 
coupling. 

The order parameter in Eq. (|40|, however, supports addi- 
tional stable topological defects .^^^^ "^^ This is tied to the fact 
that i^tripiet is invariant under combined shifts of (/) ^ (/) + tt 
and d — d, which allows for ^ half quantum vortices in 
which the superconducting phase (j) and d both rotate by tt 
around a vortex core. As a concrete example, consider the 
order parameter configuration^^ 







-ie/2 



d(r) = cos((9/2)x + sin((9/2)y, (41) 



where (r, 6) are polar coordinates. Inserting this form into Eq. 
( |40| ), one finds 

d'^r—[tlj^{dx ^idy)'il)^ 

- e-''ij^{d,^idy)^^^H.c.l (42) 

revealing a key feature of half quantum vortices — these de- 
fects are equivalent to configurations in which only one spin 
component 'sees' an ordinary |^ vortex. Thus a half quan- 
tum vortex binds a single zero-energy Majorana mode, just as 
for vortices in the spinless p -\- ip superconductor discussed 
earlier. Typically, however, nucleating half quantum vortices 
costs more energy than ordinary |^ vortices due to spin-orbit 
coupling, though clever routes of avoiding this outcome have 
been proposecP^lCll. in fact evidence of half quantum vortices 
in mesoscopic Sr2Ru04 s ample s was very recently reported 
experimentall>^^ (see Sec. IV C). 



Finally, we note in passing that it is also in principle pos- 
sible for a time-reversal-invariant 2D superconductor to form 
such that one spin undergoes p -Mppairing while its Kramer's 
partner exhibits p — ip pairin^^JES Provided time-reversal 
symmetry is present, such phases support stable counter- 
propagating chiral Majorana modes at the boundaries between 
topological and trivial regions. These can be viewed as a 
superconducting analog of 2D topological insulators, where 
counter-propagating edge states formed by Kramer's pairs are 
similarly stable due to time-reversal symmetry^. 



III. PRACTICAL REALIZATIONS OF MAJORANA 
MODES IN ID p-WAVE SUPERCONDUCTORS 

A. Preliminary Remarks 

We will now survey several ingenious schemes that have 
been proposed to realize Majorana fermions in topological 
phases similar to that of Kitaev's model for a ID spinless p- 



wave superconductor reviewed in Sec. |II A[ To put the prob- 
lem in perspective, it is useful to highlight the basic chal- 
lenges involved in realizing Kitaev's model experimentally. 
First, there is a 'fermion doubling problem' of sorts that must 
be overcome — since electrons carry spin- 1/2 one must freeze 
out half of the degrees of freedom so that the ID system ap- 
pears effectively 'spinless'. Stabilizing p- wave superconduc- 
tivity for such a 'spinless' system poses a still more serious 
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challenge. Not only are p-wave superconductors exceedingly 
rare in nature, but an attractively interacting ID electron sys- 
tem that conserves particle number can at best exhibit power- 
law superconducting correlations in contrast to the long-range 
ordered superconductivity assumed in Kitaev's model. (Re- 
markably, power-law superconducting order can be sufficient 
to stabilize Majorana modes^^ though the splitting of the 
degenerate ground states in such cases scales as a power-law 
of the system size rather than exponentially.) The proposals 
we review below employ the same three basic ingredients to 
cleverly overcome these challenges: superconducting proxim- 
ity effects, time-reversal symmetry breaking, and spin-orbit 
coupling. 

The essence of the first ingredient is that a ID system can 
inherit Cooper pairing from a nearby long-range-ordered su- 
perconductor. Fluctuations of the resulting superconducting 
order parameter for the ID system are largely controlled by 
the parent bulk superconductor, and can thus remain unimpor- 
tant even at finite temperature despite the low dimensionality 
of the parasitic material. Since superconducting proximity ef- 
fects are central to much of this review, we will digress briefly 
to elaborate on the physics in greater detail. Consider for the 
moment some ID electron system with a Hamiltonian of the 
form 



Hid 



dk 
2^ 



(43) 



and a conventional bulk s-wave superconductor described by 

Here ?/^^^ and r^^j^ add electrons with spin a to the ID 
system and superconductor, respectively, while esc{k) = 
k'^/{2msc) — l-^sc and A^c are the superconductor's kinetic en- 
ergy and pairing amplitude. When the ID system is brought 
into intimate contact with the superconductor [as in Fig.[6|a)], 
the resulting structure can be described by 



H = Hid + Hsc + ^r, 



(45) 



where Hy encodes single-electron tunneling between the two 
subsystems with amplitude F. Taking the ID system to lie 
along the line {x^y^z) = (x, 0, 0), one can explicitly write 



(46) 



The effect of the hybridization term Hr can be crudely de- 
duced using perturbative arguments and dimensional analy- 
sis. Suppose that the superconductor's Fermi wavevector 
greatly exceeds that of the ID system. Intuitively, in this 
regime (which is relevant for all of the setups of interest) the 
hybridization between the two subsystems should be primarily 
controlled by F and properties of the superconductor. When 
F/c|? <C Age, it suffices to treat Hr perturbatively since in 
this limit single electron tunneling is strongly suppressed due 
to the parent superconductor's gap. At second order one gen- 



erates an effective Cooper-pair hopping term which, using di- 
mensional analysis, takes the form 



SH ( 



hsc/\ 



j dx (V'T-V';x^I(^,o,o)^|(^,o,o) + -H"-c-) • 



At low energies one can replace 77|?7| — ?• {i^jjil^) oc p. 



(47) 

where the brackets denote a ground state expectation value 
and psc is the superconductor's density of states at the Fermi 
level. In this way one arrives at the following effective Hamil- 
tonian for the ID system. 



= Hid + 
Ha = A J dx {tp^xi^ix + H.c.) . 



(48) 



with A (X 



k-/Az(psc^sc) OC p2D^ andp2D = msc/i^n) 
the superconductor's 2D density of states atk^ = 0. 

The treatment above captures a simple effective Hamilto- 
nian for the ID system that incorporates proximity-induced 
pairing. Similar models appear frequently in the literature and 
will be employed often here as well. Several authors have, 
however, emphasized the need to treat the proximity effect 
more rigorously to obtain a quantitative understanding of the 
devices we will explore below^^ A more accurate way for- 
ward involves constructing the Euclidean action correspond- 



ing to H in Eq. ( 45 ) and then integrating out the parent su- 
perconductor's degrees of freedom. Appendix [A| sketches the 
calculation and yields the following effective action for the ID 
system. 



eff 



/duj dk 

A5c[l - Z(cj)][V^^(fc,^)V^;(_/, 



iuj + Z(a;)H/c]V^(/c,u;) 

^H.c.]} (49) 



As in our perturbative analysis, an effective Cooper pairing 
term (now frequency dependent) once again appears. This 
more rigorous procedure, however, reveals that the tunneling 
F also generates a reduced quasiparticle weight Z{uj) for elec- 
trons in the ID system given approximately b}I^ 



AL 



(50) 



The physics underlying Eqs. ( [49] ) and ( [5Q| ) is that by enhanc- 
ing F the wavefunctions for electrons in the ID system bleed 
farther into the parent superconductor, thereby reducing their 
quasiparticle weight Z{uj) and enhancing the pairing ampli- 
tude that they inherit [which can reach a maximum of A^c as 
Z{uj) 0]. The reduced quasiparticle weight also, however, 
effectively rescales the original Hamiltonian Hk and dimin- 
ishes the energy scales intrinsic to the ID system.!^ [Poles 
in the electron Green's function follow from Z(uj)l-Lk, rather 
than H/e.] In other words, in an effective ID description of the 
hybrid structure, parameters such as spin-orbit coupling. Zee- 
man splitting, etc. do not take on the values one would mea- 
sure in the absence of the superconductor, but rather are renor- 
malized downward due to the hybridization. This aspect of 
the proximity effect is often neglected, but as we will see later 
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can lead to important and counterintuitive consequences. We 
should note that even at this level the modeling of the prox- 
imity effect remains rather crude. More sophisticated treat- 
ments where one treats the pairing self-consistently are also 
possible^^ but will not be discussed here. 

Remarkably, most proposals for engineering Kitaev's 
model for a ID spinless p-wave superconductor in fact ex- 
ploit proximity effects with ordinary 5-wave superconductors 
like we treated above. (It is hard to overemphasize the im- 
portance of this feature insofar as experimental prospects are 
concerned, given the many thousands of known 5-wave super- 
conductors.) While this naively appears somewhat paradoxi- 
cal, spin-orbit coupling — typically in conjuction with time- 
reversal symmetry breaking — can effectively convert such a 
ID system into a p-wave superconductor. We will now ex- 
plore a variety of settings in which such a mechanism appears. 



B. 2D Topological Insulators 

In 2005, a revolution in our understanding of a seemingly 
well-understood phase of matter — the band insulator — began 
to emerge^ ^ It is now appreciated that such states need 

not be trivial in the sense of having no available low-energy 
degrees of freedom at zero temperature. Rather, there ex- 
ists a class of topological band insulators that while inert in 
the bulk necessarily possess novel conducting states at their 
boundary. These topological phases can appear in either two- 
or three-dimensional crystals and, remarkably, merely require 
appreciable spin-orbit coupling and time-reversal symmetry. 
The numerous fascinating developments that grew out of the 
discovery of topological insulators include Fu and Kane's pi- 
oneering proposals^^ for generating Majorana fermions at 
their edges (in 2D crystals) or surfaces (in 3D). In this section 
we will describe how one can engineer a topological super- 
conducting state similar to that of Kitaev's model using the 
edge of a 2D topological insulator; the 3D case will be re- 



viewed in Sec. |IVD[ (See also Sec. [Ill D| for a proposal in- 
volving nano wires built from 3D topological insulators.) 

The hallmark of a 2D topological insulator is the presence 
of counter-propagating, spin-filtered edge states that are con- 
nected by time-reversal symmetry. In an oversimplified pic- 
ture that is adequate for our purposes, one can envision spin up 
electrons propagating clockwise around the edge while their 
Kramer's partners with spin down circulate counterclockwise 
as shown in Fig.jSja). These low-energy edge modes can be 
described by the Hamiltonian 



2DTI 



(51) 



where v is the edge-state velocity, fi is the chemical potential, 
and ijj'l^ adds an electron with spin cr at position x along the 
edge. The blue and red lines of Fig. |5jb) sketch their disper- 
sion. Provided time-reversal symmetry is preserved (elastic) 
backscattering between the counter-propagating edge states is 
prohibited even in the presence of strong non-magnetic disor- 
der. Consequently these modes are robust against localization 
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FIG. 5. (a) Schematic of counter-propagating, spin-filtered edge 
states in a 2D topological insulator, (b) Edge- state dispersion when 
time-reversal symmetry is present (red and blue lines) and with a 
Zeeman field h of the form in Eq. ([55]) (solid curves), (c) A proxi- 
mate s-wave superconductor drives the edge into a topological phase 
similar to the weak-pairing phase in Kitaev's toy model for a ID 
spinless p-wave superconductor. When the Zeeman field h is present, 
the topological phase survives provided h < y^A^ + /i^ leading to 
the phase diagram in (d). Domain walls between topological (green 
lines) and trivial regions (dashed lines) on the edge trap localized 
Majorana zero-modes. As described in the text these can be created 
with (e) a ferromagnetic insulator, (f) a Zeeman field combined with 
electrostatic gating, or (g) applying supercurrents near the edge. 



that plagues conventional ID systems. Furthermore, by fo- 
cusing on these edge states one immediately beats the fermion 
doubling problem noted earlier — the spectrum supports only a 
single pair of Fermi points as long as the Fermi level does not 
intersect the bulk bands, and in this sense the system appears 
'spinless'. 

Realizing topological superconductivity then simply re- 
quires gapping out the edge via Cooper pairing. Since the 
counter-propagating edge modes carry opposite spins, this can 
be achieved by interfacing the topological insulator with an or- 
dinary 5-wave superconductor^^; see Fig.jsj^c). As discussed 
in Sec. [Ill A| the superconducting proximity effect on the edge 
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can be crudely modeled with a Hamiltonian 

H = i^2DTI + (52) 
Ha = j dxAi^^^^^H.c), (53) 

where A is the pairing amplitude inherited from the nearby 
superconductor. Equation ([52]) yields quasiparticle energies 

(54) 



with k the momentum, and describes a gapped topological 
superconductor that is a time-reversal-symmetric relative of 
the weak-pairing phase in Kitaev's toy model. Let us now 
make this connection more precise and elucidate how the spin- 
singlet pairing A mediates p-wave superconductivity. 

To this end it is instructive to violate time-reversal symme- 
try by introducing a Zeeman field that cants the spins away 
from the z direction: 



Hz 



(55) 
(56) 



^2DTI + Hz + H/\ 

-h J dxij)^ a'' 

with h > {) the Zeeman energy. When A = the edge- 
state spectrum becomes e±(k) = —ji ± ^/^(vky^^^^}? , and 
as shown by the solid black lines in Fig. [^b) exhibits a gap 
at /c = due to the broken time-reversal symmetry. To un- 
derstand the influence of proximity-induced pairing, we first 
note that the effect of A is obscured by the fact that Eq. ( [55] ) 
contains a standard spin- singlet pairing term but an uncon- 
ventional kinetic energy form (due to the interplay of spin- 
momentum locking and the field). The physics becomes much 
more transparent upon expressing H' in terms of operators 
tl^l^{k) that add electrons with energy e±{k) to the edge. In 
this basis H' reads 



I 



dk 
27r 



{e+{k)i,X{k)i,+ {k) + e_(fc)Vl(fc)V-(fc) 
[',p+{-k)^P+{k) + i^-{-k)^P-{k) + H.c] 



+ A,(fc)[^/;_(-fc)V+(fc) + J?.c.]}, 
where the pairing functions are 



vkA 



A,(fc) 



hA 



(57) 



(58) 



The first line of Eq. ( [57] ) simply describes the band energies 
while the third captures interband 5-wave pairing. Most im- 
portantly, the second line encodes intraband p-wave pairing. 
This emerges because, as shown schematically in Fig. [5jb), 
electrons at k and — /c in a given band have misaligned spins 
and can thus form Cooper pairs in response to A. By Fermi 
statistics, the effective potential Ap{k) that pairs these elec- 
trons must exhibit odd parity since they derive from the same 
band. (Physically, the odd parity reflects the fact that the elec- 
tron spins rotate as one sweeps the momentum from kto —k.) 
This is the first of many instances we will encounter in which 
an 5-wave order parameter effectively generates p-wave pair- 
ing by virtue of spin-orbit coupling. 



The connection to Kitaev's model becomes explicit in the 
limit where h A and /i resides near the bottom of the upper 
band as in Fig.[5jb). In this case the lower band plays essen- 
tially no role and can be projected away by simply sending 
?/^_ ^ 0. Furthermore, only momenta near k = are impor- 
tant here so it suffices to expand e+ (/c) ^ —{ii — h)^^k'^ = 
-/ieff + A:^/(2meff) and Ap{k) ^ ^k = A^i^k. With these 
approximations, one obtains an effective Hamiltonian that in 
real space reads 



^eff 



dx 

Aeff 



2meff 



Meff 



(59) 



which describes Kitaev's model for a ID spinless p-wave su- 
perconductor in the low-density limit \i.e., near /i = —I'm 
Fig.[TJa)]. A similar mapping can be implemented for /i near 
the top of the lower band in Fig. [5jb). These considerations 
show that for ^ A, the edge forms a trivial strong pairing 
phase when \pi\^h and a topological weak pairing phase at 

ImI > h. 

A more accurate phase diagram valid at any /i, A can be 



deduced by studying the unprojected Hamiltonian in Eq. ( 57 ), 
which yields quasiparticle energies 



E'^(k) 



A2. 



±(e+-e_)yA2T^.(60) 



The quasiparticle gap extracted from Eq. ( 60 ) vanishes only 
when = A^ + /i^. Matching onto the A results 

derived above, we then conclude that the edge forms a topo- 
logical superconductor provided 



h < \/A^ + /i^ (topological criterion). (61) 

Physically, the edge forms a topological phase if supercon- 
ductivity dominates the gap but a trivial phase if the gap is 
driven by time-reversal symmetry breaking. Figure [5jd) illus- 
trates the resulting phase diagram. Note that topological su- 
perconductivity persists even in the time-reversal- symmetric 
limit with h = this has important physical consequences 
as we discuss shortly. We also note that electrons on the 
edge are additionally subject to Coulomb repulsion, which 
have dramatic consequences in ID and have so far been ne- 
glected. References [99] and 1 1001 find that while strong inter- 
actions (with a Luttinger parameter g < 1/2) can destroy the 
topological phase, milder repulsion (1/2 < g < 1) leaves the 
phase diagram of Fig.[5jd) qualitatively intact. 

Thus far we have only shown how to utilize the edge to 
construct a ID topological superconductor on a ring, with no 
ends. Stabilizing localized Majorana zero-modes requires in- 
troducing a domain wall between gapped topological and triv- 
ial phases on the edge^^. The setups of Figs.[5];e) and (f) trap 
Majorana modes 71,2 by gapping three sides with supercon- 
ductivity and the fourth with a Zeeman field h of the form 



in Eq. ( [56] ). Topological and trivial regions are respectively 
indicated by green and dashed lines in these figures. In (e) 
electrons on the lower edge 'inherit' the Zeeman field via a 
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proximity effect with a ferromagnetic insulator, just as a pair- 
ing field A is inherited from a superconductor.^^ Note that the 
chemical potential for the bottom edge must reside within the 
field-induced spectral gap [recall Fig.[5jb)]; otherwise that re- 
gion remains gapless despite the broken time-reversal. In (f) 
both superconductivity and the Zeeman field h are uniformly 
generated on all four edges, the latter by applying a magnetic 
field. Provided h > A, the topological and trivial regions 
form simply by adjusting the chemical potential /i via gating 
so that h < Y^A^ + /i^ on three sides while h > y^A^ + /i^ 
on the fourth. We emphasize here that one can simultane- 
ously have h > A and still be well below the superconduc- 
tor's critical field since the Zeeman energy for the edge can 
greatly exceed that in the superconductor due to spin-orbit en- 
hancement of the ^-factoi^^^. 

Majorana zero-modes can also be trapped by selectively 
driving supercurrents near the edge of the sample To un- 
derstand the underlying principle, let us revisit the Hamilto- 
nian in Eq. ( [55] ) when a supercurrent / flows as in Fig. [5jg). 
The current generates a phase twist in the superconducting or- 
der parameter so that H/\ becomes 



Ha 



(62) 



with / (X dxcj^ix). It is convenient to gauge away the phase 
factor above by sending tpcr e~*'^'^^^/^?/^cr; defining hz = 
vdx(l)/2, the full Hamiltonian then reads 



dx{ilj^ [-{ivdx + hz)a^ - M - ha^] i/j 



(63) 



Equation ( [63] ) shows that the supercurrent mimics the effect 
of a Zeeman field hz directed along the z direction. Con- 
trary to h, this does not open a gap but rather breaks the res- 
onance between electrons with momentum k and —k, thereby 
suppressing their ability to Cooper pair. Suppose now that 
< h < V'A^ + When / = the edge then forms 
a topological phase where A dominates the gap. At large / 
(such that hz ^ A), however, the pair-breaking effect of hz 
essentially kills A and the edge forms a trivial state with a 
gap arising from h. Supercurrents therefore allow one to turn 
a topological portion of the edge into a trivial state, similar to 
the gate in Fig.[5jf), providing yet another means for stabiliz- 
ing Majorana-carrying domain walls. 

In our view 2D topological insulators hold great promise 
as a potential venue for Majorana fermions, particularly in 
the long term. For one, their reduced dimensionality should 
allow for bulk carriers — which usually bedevil 3D topolog- 
ical insulators — to be removed relatively easily by electro- 
static gating. The topological superconducting phase hosted 
by the edge also exhibits several remarkable features. First, 
this phase is 'easy' to access in the sense that its appearance 
requires the chemical potential /i to satisfy the inequality in 



Eq. (61) while not intersecting the bulk bands; this chemi- 
cal potential window is therefore limited by the bulk gap for 
the topological insulator which can reach the ~ O.leV scale 
(see, e.g., Ref. 1 1041) . By contrast the trivial gapped state re- 
quires positioning /i inside of the Zeeman-induced gap of Fig. 



[5jb), which likely requires greater care. While ultimately the 
ability to access both kinds of states is essential, the compar- 
ative ease for forming the topological phase greatly facilitates 
the Josephson-based Majorana detection schemes discussed 
in Sec. [VB[ More strikingly, as a consequence of Ander- 
son's theorem the gap protecting the time-reversal-invariant 
topological superconductor that forms when h = is unaf- 
fected by non-magnetic disordeiElMI, We emphasize that one 
needn't work at = to enjoy this protection: with h ^ but 
/i far from the Zeeman-induced gap, electrons near the Fermi 
energy are weakly perturbed by the field and hence 'almost' 
obey Anderson's theorem^"^. 

A final noteworthy feature pertains to how large the topo- 
logical superconductor's gap can be in principle. Addressing 
this question requires the more rigorous treatment of the prox- 



imity effect discussed in Sec. Ill A Recall that increasing the 
tunneling F between the superconductor and topological insu- 
lator enhances A but reduces the energy scales intrinsic to the 
edge. When h = [such as in the setup of Fig. [5je)] it fol- 
lows from Eq. ( 54 ) that the gap is simply £^gap = A, which is 
independent of the quantities /i that F suppresses. Thus in 
this case the gap increases monotonically with F, reaching a 
maximum of A^c for the parent superconductor.'^^'^ In other 
words, it is in principle possible for the edge to inherit the full 
pairing gap exhibited by the parent superconductor. This is- 
sue becomes subtler in setups such as Fig.[5jf) where h ^ 0, 
for in this case increasing F supresses the Zeeman-induced 
gap in the spectrum, making it more difficult to stabilize the 
trivial phase to trap Majoranas. How large a hybridization is 
desirable then depends on details such as the tolerable fields 
one can apply, sample purity, etc. 

Despite these virtues this platform for Majorana fermions 
faces the hurdle that experimental progress on 2D topolog- 
ical insulators has to date remained rather limit ed. Though 
numerous candidate materials have been proposecP^E^^ISltnni 
only predictions for HgTe have so far been confirmed 
experimentall}Wii^. Some evidence for a topological in- 
sulator phase in InAs/GaSb quantum wells also appeared 
recentljC^^C^, though the signatures are less clear cut due 
to persistence of bulk carriers in the samples. The situa- 
tion, however, already shows signs of improvement — very re- 
cently topological insulator behavior in HgTe has been inde- 
pendently confirmed by the Yacoby group, and experimental 
efforts to introduce a proximity effect at the edge are under- 
way. It will be very interesting to see how this avenue pro- 
gresses in the near future. 



C. Conventional ID wires 



Two seminal works (Lutchyn et al^^ mid Oreg et al!^^^) re- 
cently established that one can engineer the topological phase 
in Kitaev's toy model by judiciously combining three exceed- 
ingly simple and widely available ingredients: a ID wire with 
appreciable spin-orbit coupling, a conventional 5-wave super- 
conductor, and a modest magnetic field. Figure[6|;a) illustrates 
the basic architecture required, which can be modeled by the 
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s-wave superconductor 



FIG. 6. (a) Basic architecture required to stabilize a topological superconducting state in a ID spin-orbit-coupled wire, (b) Band structure 
for the wire when time-reversal symmetry is present (red and blue curves) and broken by a magnetic field (black curves). When the chemical 
potential lies within the field-induced gap at /c = 0, the wire appears 'spinless'. Incorporating the pairing induced by the proximate super- 
conductor leads to the phase diagram in (c). The endpoints of topological (green) segments of the wire host localized, zero-energy Majorana 
modes as shown in (d). 



following Hamiltonian, 



H — i^wire + H/\ 



(64) 



/ ^"^^^ l~S ^o^^^dx + ha' ) V^(65) 



A = y (ixA(V^^V^; + H.C.). 



(66) 



Here adds an electron with effective mass m, chemical po- 
tential /i, and spin a to the wire; a > denotes the strength for 
spin-orbit coupling that favors aligning spins along or against 
the y direction depending on the momentum; and > is the 
Zeeman energy arising from a magnetic field applied along 
z. (The precise spin-orbit and magnetic field axes are unim- 
portant so long as they are perpendicular.) For concreteness 
one can envision i^wire describing an electron-doped semi- 
conducting wire such as In As with Rashba coupling in 
the limit where only the lowest transverse subband is rele- 
vant. The pairing term H/\ crudely models the proximity ef- 
fect on the wire arising from the adjacent 5- wave supercon- 
ductor. Note that H above takes on the same form as the 



topological insulator edge Hamiltonian in Eq. (55), with 
the sole addition of an ordinary k'^ / (2m) kinetic energy con- 
tribution. This important modification underlies many of the 
qualitative distinctions between the two setups. 

Let us first consider A = and elucidate how the Hamilto- 
nian i^wire overcomes our fermion doubling problem. The red 
and blue curves in Fig. [6jb) illustrate the wire's band struc- 
ture in the limit where h = 0. Due to spin-orbit coupling, 
the blue and red parabolas respectively correspond to elec- 
tronic states whose spin aligns along -\-y and —y. Clearly no 
'spinless' regime is possible here — the spectrum always sup- 
ports an even number of pairs of Fermi points for any /i. The 
magnetic field remedies this problem by lifting the crossing 
between these parabolas at A: = 0, producing band energies 



2m 



(67) 



sketched by the solid black curves of Fig. [6jb). When the 
Fermi level resides within this field-induced gap (e.g., for /i 
shown in the figure) the wire appears 'spinless' as desired. 

The influence of the superconducting proximity effect on 
this band structure can be intuitively understood by focusing 



on this 'spinless' regime and projecting away the upper unoc- 
cupied band, which is legitimate provided A <C /i. Crucially, 
because of competition from spin-orbit coupling the mag- 
netic field only partially polarizes electrons in the remaining 
lower band as Fig.[6jb) indicates schematically. Turning on A 
weakly compared to h then effectively p-wave pairs these car- 
riers, driving the wire into a topological superconducting state 
that connects smoothly to the weak-pairing phase of Kitaev's 
toy model (see Ref.[35lfor an explicit mapping). 

More formally, one can proceed as we did for the topolog- 
ical insulator edge and express the full, unprojected Hamil- 
tonian in terms of operators V^±(/c) that add electrons with 
energy e±{k) to the wire. The resulting Hamiltonian is again 
given by Eqs. ( [57]) an d ( 58 ) [but with v ^ a and band energies 
e±{k) from Eq. (|67|)], explicitly demonstrating the intraband 



p-wave pairing mediated by A. Furthermore, Eq. (60) pro- 
vides the quasiparticle energies for the wire with proximity- 
induced pairing and again yields a gap that vanishes only 
when h = y^A^ + /i^. For fields below this critical value the 
wire no longer appears 'spinless', resulting in a trivial state, 
while the topological phase emerges at higher fields. 



h> \/A2 + /i^ (topological criterion). 



(68) 



Figure [6jc) summarizes the phase diagram for the wire. No- 
tice that this is inverted compared to the topological insulator 
edge phase diagram in Fig. [5jd). This important distinction 
arises because the k'^/ (2m) kinetic energy for the wire causes 
an upturn in the lower band of Fig.[6];b) at large thereby 
either adding or removing one pair of Fermi points relative to 
the edge band structure. 

Since a wire in its topological phase naturally forms a 
boundary with a trivial state (the vacuum), Majorana modes 
7i and 72 localize at the wire's ends when the inequality in Eq. 



(68) holds. Majorana- trapping domain walls between topo- 
logical and trivial regions can also form at the wire's interior 
by applying gate voltages to spatially modulate the chemical 
potential^^ or by driving supercurrents through the adja- 
cent superconductor^^^ (using the same mechanism discussed 
in Sec. |IIIB| ). Figure [6];d) illustrates an example where four 
Majoranas form due to a trivial region in the center of a wire. 

It is useful address how one optimizes the ID wire setup 
to streamline the route to experimental realization of this pro- 
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posal. This issue is subtle, counterintuitive, and difficult even 
to define precisely given several competing factors. First, how 
well should the wire hybridize with the parent superconduc- 
tor? The naive guess that the hybridization should ideally be 
as large as theoretically possible to maximize the pairing am- 
plitude A imparted to the wire is incorrect. One practical issue 
is that exceedingly good contact between the two subsystems 
may lead to an enormous influx of electrons from the super- 
conductor into the wire, pushing the Fermi level far above the 
Zeeman-induced gap of Fig.[6jb) where the topological phase 
arises. Restoring the Fermi level to the desired position by 
gating will then be complicated by strong screening from the 
superconductor. 

Reference [94l emphasized a more fundamental issue related 
to the optimal hybridization. The topological phase's sta- 
bility is determined not only by the pairing gap induced at 
the Fermi momentum, Ek^^ cx A, but also the field-induced 
gap at zero momentum, Eq = \h — y^A^ + required to 
open a 'spinless' regime. The minimum excitation gap for 
the topological phase is set by the smaller of these two en- 
ergies. As reviewed in Sec. |III A| increasing the tunneling F 
between the wire and superconductor indeed enhances A but 
simultaneously reduces the Zeeman energy h. From the ef- 
fective action in Eq. ( 49 ) we explicitly have h = Z/ibare and 
A = (1 — Z)Asc, where /ibare is the Zeeman energy for the 
wire when the superconductor is absent, A^c is the parent su- 
perconductor's gap, and Z is the quasiparticle weight defined 
in Eq. ( [50| (for simplicity we neglect the frequency depen- 
dence). Suppose now that F increases from zero, thereby re- 
ducing Z from unity. As Z decreases the topological phase's 
gap initially increases due to an enhancement of Ekp . Even- 
tually, however, the gap decreases due to a suppression of Eq, 
and the topological phase disappears entirely beyond a critical 
value of F at which Eq vanishes. The maximum achievable 
gap depends sensitively on details such as the spin-orbit cou- 
pling strength, applied field amplitude, mobility for the wire, 
etc^ 

But should one ideally design the setup to achieve this max- 
imum gap? This, too, is not necessarily the case. As an il- 
lustrative example, consider a 200, OOOcm^A^s mobility wire 
with parameters relevant for InAs, adjacent to a superconduc- 
tor with Age = 2K. Reference l97l predicted (including short- 
range disorder) that such a wire realizes an optimized gap 
of £^max ^ 0.3K when the Zeeman energy is /ibare ^ IK. 
Realizing the topological phase in a meaningful way then re- 
quires positioning the chemical potential within a rather nar- 
row ~ IK window over distances long compared to the wire's 
coherence length, which could prove challenging experimen- 
tally due to disorder-induced chemical potential fluctuations. 
Thus it may be desirable to apply larger magnetic fields to 
soften these constraints at the expense of reducing the gap 
somewhat. 

The question of how large this Zeeman field should be is 
also rather delicate. On one hand, enhancing h indeed alle- 
viates the need to fine tune fi. But on the other, increasing 
h suppresses superconductivity in the parent superconductor, 
reduces the effective p-wave pairing amplitude for electrons 
in the wire due to a further alignment of their spins [see Eq. 



(58)], and makes the wire more susceptible to disordei^^^H. 
As discussed in the topological insulator context, the first ef- 
fect can be rather minor if the g factor in the wire greatly 
exceeds that in the superconductor. The effect of disorder 
warrants more serious consideration. Since the topological 
phase appears only at finite magnetic fields Anderson's the- 
orem does not protect the gap against disorder in the wire — 
which is al ways pai r-breaking in this context as many studies 
have showiJS^^^. (Fortunately though, Refs. [I26| and 
"127 'conclude that disorder native to the proximate 5- wave su- 
perconductor is benign.) 

The sensitivity of the topological phase to disorder is de- 
termined by how severely time-reversal symmetry is bro- 
ken. One can quantify this by the ratio of the Zeeman en- 
ergy h to the spin-orbit energy which we define here by 
Eso = \ma^. Physically, Ego is the Fermi energy measured 
relative to the bottom of the bands when /i = = 0; see 
Fig.[6jb). At large h/Eso spins near the Fermi level are fairly 
well polarized, so time reversal is strongly violated and hence 
disorder can efficiently suppress the gap. In contrast, at small 
h/Eso spins at kp and —kp are nearly antiparallel due to the 
dominance of spin-orbit coupling. Time-reversal is then 'al- 
most' present insofar as carriers near the Fermi level are con- 
cerned, thereby sharp ly sup pressing the impact of disorder on 
the topological phase. 1^^^ 

The strength of spin-orbit coupling is thus a crucial materi- 
als parameter. 'Large' spin-orbit values allow one to operate 
at relatively high Zeeman fields — where the topological phase 
occurs over a broad chemical potential range — while main- 
taining some degree of robustness against disorder. (Though 
one should bear in mind that increasing spin-orbit coupling 
can reduce the mobility, thus at least partially offsetting this 
advantage. ^^) Furthermore, the maximum gap that the topo- 
logical phase can in principle exhibit increases with the spin- 
orbit strength,^"^ approaching a value of A^c in the limit 
where^"^ E^o > > A^c- 

In light of this discussion it is interesting to ask how 
electron-doped InAs and InSb wires fare as platforms for Ma- 
jorana fermions. At present these are the most commonly dis- 
cussed wires for this proposal, and for good reason. Both ex- 
hibit exceptionally large g factors (^inAs ^15 and ^inSb ^ 50 
for bulk crystals), and can be synthesized with high mobility 
and long mean free paths. Good superconducting proximity 
effects have also been measured in both systems. ^^^ ^^^ One 
challenge with these materials, however, is that while they are 
often lauded as having strong spin-orbit coupling, the energy 
scale Eso is typically of order IK for both InAs and InSb^. 
Rashba coupling is gate-tunable to some extent,E^but it may 
nonetheless prove difficult to access the topological phase in 
a spin-orbit dominated regime where h/Ego <C 1. Disorder 
is thus likely to play a nontrivial role in these settings. Still, 
it is hard not to be optimistic about the prospects of success 
given the high level and rapid pace of ongoing experimental 
activity. 

Several subsequent works have pursued variations on this 
proposal in an effort to mitigate the challenges involved with 
realizing the topological phase experimentally. One issue that 
has been explored is whether one can reduce the applied mag- 
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netic field required to access the topological phase. Refer- 
ence |99|proposed the use of nuclear spins to generate the Zee- 
man energyES, thus removing the external magnetic field alto- 
gether. Repulsive interactions — which are inevitably present 
in the wire — also allow the topological phase to be accessed 
at weaker magnetic fields and over a broader chemical poten- 
tial window due partly to exchange enhancement of the Zee- 
man field^^"^ (but if the repulsion is too strong this state can 
disappea r^^P^ . In principle, interactions can even stabilize 
topological superconductivity at zero magnetic field due to 
spontaneous time-reversal symmetry breaking. ^^"^ Generating 
the topological phase by applying supercurrents as described 
in Ref. 103 also lowers the critical Zeeman field needed. 

Multichannel wires have been shown by numerous studies 
to support a ID topological superconducting state away from 
the lowest-subband limit.^^^ ^^"^ ^"^^ In this more general 
setting one simply needs an odd number of partially occupied 
bands and a wire whose width does not exceed the coherence 
length. The former criterion is rather natural: starting from the 
topological phase in a single-channel wire, pairs of partially 
occupied bands can always be adiabatically introduced with- 
out closing the gap. The latter criterion ensures that the in- 
duced superconducting phase remains quasi- ID, which is re- 
quired for the wire to exhibit a substantial gap.E^^CMnilxhese 
works are significant in part because multiple subbands are 
usually occupied in semiconducting wires such as InAs and 
InSb. Gating these wires into the lowest subband regime may 
be nontrivial particularly when a superconductor is nearby, but 
fortunately is unnecessary. Furthermore, multichannel wires 
open the door to realizing a topological phase in a variety of 
other settings, such as gate-defined channels in quantum wells 
or surface states featuring large spin-orbit coupling^^^*^^^^^. 

Even in a multichannel wire accessing the topological 
phase will require some degree of gating. An interesting pos- 
sible route to enhancing gate-tunability is to employ periodi- 
cally modulated structures in which a regular array of super- 
conducting islands contacts the wire.'^^^'SMl xhe basic idea is 
that gating the wire in the regions between adjacent super- 
conductors may be relatively easy. Although the Hamilto- 
nian is rather different from that of the uniform structure con- 
sidered previously, a topological phase can still arise, which 
can be understood simply in two limits. First, suppose that 
the chemical potential varies only modestly along a single- 
channel wire. In the limit where the Fermi wavelength ex- 
ceeds the spacing between superconducting islands the elec- 
trons effectively 'see' only the average induced pairing poten- 
tial. Periodic modulations are then essentially smeared out, 
and a topological phase arises under similar conditions to the 
uniform case. 

Second, suppose that a large potential barrier formed at the 
boundary between the gated and superconducting regions of 
the wire, effectively creating a chain of quantum dots bridged 
by superconducting islands. ^"^^ By introducing a magnetic 
field and fine-tuning the gates, one can in principle bring a sin- 
gle level into resonance on each dot. Effectively each dot then 
behaves as a single site in a 'spinless' chain. Electrons from 
neighboring sites communicate indirectly via the supercon- 
ductor, which can mediate both the nearest-neighbor hopping 



and p-wave pairing in Kitaev's toy model. This setup can be 
particularly promising if the hopping and pairing amplitudes 
can be tuned equal to one another; this limit corresponds to the 
t = A case discussed in Sec. II A where Majorana end-states 
localize at a single site. Majorana zero-modes might then be 
observable in an array consisting of relatively few quantum 
dots. One obvious challenge here is the high level of fine- 
tuning required to reach this regime. Furthermore, strong ran- 
domness may pose an issue given that the hopping and pairing 
parameters presumably depend exponentially on factors such 
as the width of the superconducting islands, barrier heights, 
etc. 

Numerous other interesting variants have been introduced. 
Hole-doped semiconducting wires — which benefit from a 
greatly enhanced spin-orbit energy Ego relative to their 
electron-doped counterparts — comprise one ver y promising 
alternative for realizing Majorana modes .^322110 Carbon nan- 
otubes can also in principle host Majoranas despite the fact 
that obtaining a 'spinless' regime with proximity-induced 
pairing is rather nontrivial fl^iiU^ Another interesting proposal 
involves a half-metallic ferromagnetic wire in which only one 
spin species conducts. While a 'spinless' regime emerges triv- 
ially here, our standard trick for inducing p-wave supercon- 
ductivity via a conventional s-wave superconductor no longer 
works. (A spin-singlet order parameter can not pair spins that 
are fully aligned). This problem can be solved by coupling 
the half-metal to a non-centrosymmetric superconductor with 
spin-orbit coupling. Such a superconductor generically con- 
tains both spin- singlet and spin-triplet Cooper pairing and 
can therefore induce a proximity effect in the wire to generate 
topological superconductivity.QSl Remarkably, clever routes to 
engineering a topological phase in systems without spin-orbit 
coupling were even devised recently. The key idea can 
be understood by rewriting our original wire Hamiltonian in 



Eq. (64 ) in terms of rotated operators i/j = e*^"^^ i/j: 



H 



2m 



heff(x) = — sin(2max)x + cos{2max)z. 



(69) 
(70) 



It follows that the Hamiltonian for a spin-orbit-coupled wire 
subjected to a uniform magnetic field is unitarily equivalent 
to that of a spin-orbit- free wire in a spatially rotating Zeeman 
field — so either system can support a topological phase. The 
required non-uniform Zeeman field in the latter setup can be 
generated using an array of magnetic nanoparticles deposited 
on a superconductor^^^ or by coupling a set of magnetic gates 
to a wird^. Finally, we note that one can engineer a ID 
topological superconducting state with cold fermionic atoms 
using optical Raman transitions (to generate effective spin- 
orbit coupling and Zeeman fields) and a proximity effect with 
a bulk molecular BEC.^^^ This route is especially tantaliz- 
ing given the recent pioneering experiments by Lin et al]^, 
where Raman lasers created a band structure for bosons simi- 
lar to Fig.|6];b). 
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(C)$ = he/ 4:6. 



FIG. 7. (a) A 3D topological insulator nanowire with a magnetic field 
applied along its axis can realize a topological superconducting state 
when in contact with an s-wave superconductor, (b)-(d) Nanowire 
band structure when the flux $ passing through its center is (b) 0, (c) 
1^, and (d) Solid and dashed curves respectively denote doubly 
degenerate and non-degenerate bands. Green shaded regions indicate 
chemical potential windows in which an odd number of channels is 
occupied, as required for generating topological superconductivity. 



D. 3D topological insulator nanowires 

Three-dimensional topological insulators, like their 2D 
counterparts, are strongly spin-orbit-coupled materials in 
which electrons are insulating in the bulk but due to topol- 
ogy form novel metallic states at their boundary.^ ^ ^^^^ (Re- 
markably, these boundary states were captured very early on 
in Refs. ll6Ql and ll61l ) In the simplest cases each surface hosts 
a single Dirac cone described by 



H: 



3DTI 



/ 



(V X a) 



(71) 



where n is the surface normal, ip^^ adds an electron with spin 
a at position r on the surface, /i is the chemical potential, and 
V is the surface state velocity. This Hamiltonian favors orient- 
ing the electron spins along the surface, but perpendicular to 
the momentum, similar to Rashba coupling in a 2D electron 
gas. 

Cook and Franz recently showed that nanowires fabricated 
from 3D topological insulators can host a topological super- 
conducting state with several advantageous features. Con- 
sider a cylindrical nanowire of radius R, whose axis ori- 
ents along the z direction as sketched in Fig. |7ja). Rewrit- 
ing Eq. ([tT]) in cylindrical coordinates and using n = r 



for the surf ace nor mal, one obtains the following nanowire 
Hamiltonian,^!!!!] 



TI wire 



/ 



dzdOip'^ 



R 



-do — 6 • ad. 



(72) 

It is convenient to remove the angular dependence in the dz 
term by defining a new field = e*^^^/^?/^ that exhibits anti- 
periodic boundary conditions, ?/^(6> + 27r, z) = —ijj{0^z). 



After absorbing a constant into the definition of /i the Hamil- 
tonian then reads 



Hr 



TI wire 



dzdOtjj^ 



-IV I 



^ (73) 



and can be trivially diagonalized. Equation ( [73] ) admits band 
energies en±(k) = :Lv^/JrilRf^^!n^ — ji, where k is the mo- 
mentum along the cylinder axis and n is a half-integer angular 
momentum quantum number due to the boundary conditions 
on 'ij). Figure[7[b) illustrates the spectrum; each band is doubly 
degenerate so our fermion doubling problem remains here. 

Applying a magnetic field along the cylinder axis as in Fig. 
[Tj^a) heals this problem in an interesting way.^^^ Upon incor- 
porating the field by sending —ide —ide — ^/^o, where 
<l>o = hc/e and $ is the flux piercing the cylinder, the band 
energies become 



en±{k) ^ ±vy^[{n-^/^o)/R?^k^ - (74) 

Note that we have neglected the Zeeman term ^g/j^sBa^ 
since this contribution merely renormalizes ^. Generally, 
the flux lifts the band degeneracies present at zero field as 
Figs. [Tj^c) and (d) respectively illustrate for $ = |^ and 
1^ ; here dashed curves represent non-degenerate bands while 
solid curves are doubly degenerate. This produces chemical 
potential windows (shaded green regions in Fig. [7]) in which 
the electrons partially occupy an odd number of bands as de- 
sired. In these odd-channel regimes Cooper pairing states at 
the Fermi level can drive the nanowire into a ID topological 
superconducting phase. Because of spin-orbit coupling this 
can be achieved in the standard way using the proximity effect 
with an 5-wave superconductor as in Fig.jTJa). 

The case of |^ flux is particularly appealing. Here an odd 
number of surface-state bands is occupied for any ji that re- 
sides in the bulk band gap. Accessing the topological su- 
perconducting phase is then 'easy' and automatically results 
in localization of Majorana zero-modes at the ends of the 
nanowire. Recall that topological superconductivity is simi- 
larly 'easy' to obtain in a 2D topological insulator edge, but 
there the formation of Majoranas is less trivial; see Figs.[5je) 
and (f). Moreover, in the (fictitious) limit where the |^ flux 
is confined to a thin solenoid passing through the nanowire 's 
center, the bulk retains time-reversal symmetry and thus main- 
tains immunity against non-magnetic disorder. ^^^^^^ (Time- 
reversal is always lifted at the nanowire ends where the 
solenoid enters and exits the wire, which is essential for the 
formation of Majorana end- states.) In an actual experiment 
additional ingredients — such as the Zeeman effect, fluctua- 
tions in the nanowire radius that cause the flux to deviate lo- 
cally from 1^, etc. — will inevitably remove this exact protec- 
tion, though some robustness against disorder likely survives. 
Other flux values are also interesting despite the absence of 
these features. For example, near <l> = the topological and 
trivial superconducting phases appear over chemical poten- 
tial windows of roughly equal size, so tuning between these 
phases by gating may be ideal here. 

To realize this proposal the nanowire's radius R should sur- 
pass the surface state penetration depth (to have well-defined 



18 



surface states) but not exceed the induced superconducting co- 
herence length in the wire (so that the superconducting state 
remains quasi- ID). Since penetration depths of a few unit cells 
can arise ^^"^ while the coherence length is typically on the mi- 
cron scale, these criteria rather loosely constrain the wire size. 
Another relevant scale is the spacing v/R between adjacent 
bands in Fig. |7];b)-(d). Enhancing v/R allows one to sup- 
press unwanted interband coupling generated, e.g., by disor- 
der, so reducing R significantly below the coherence length 
will likely prove worthwhile. Wires with of a few tens 
of nanometers should allow one to achieve a sizable inter- 
band spacing (~ lOOK) and reach |^ flux with sub-Tesla 
fields. As experiments on topological insulator nanowires and 
nanoribbons are now steadily progressing, this approach cer- 
tainly warrants further attention. The prospect of using a weak 
magnetic field to stabilize a topological phase and Majorana 
modes without requiring careful control over the chemical po- 
tential is well worth pursuing. 



IV. PRACTICAL REALIZATIONS OF MAJORANA 
MODES IN 2D p + ip SUPERCONDUCTORS 



A. Preliminary remarks 

The proposals we review below for experimentally realiz- 
ing the physics of topological 2D p -\- ip superconductivity 
discussed in Sec. |IIB| loosely fall into two categories. The 
first corresponds to 'intrinsic' realizations, where p -\- ip pair- 
ing emerges by virtue of a material's internal dynamics. We 
briefly discuss two classic systems of this type: the fractional 
quantum Hall state at filling factor u = b/2 and the layered 
spinful triplet superconductor Sr2Ru04. 

Most of our discussion will center around the second 
category — 'engineered' topological phases — wherein 'spin- 
less' p+zp superconductivity is driven by forming appropriate 
heterostructures with various kinds of 2D electron systems. In 
this approach one faces the same basic hurdles that arose in 
realizing Kitaev's ID model: removing the spin degeneracy 
so that the system appears effectively 'spinless' and inducing 
long-range-ordered p-\-ip pairing in the remaining Fermi sur- 
face. As we will see, time-reversal- symmetry breaking, spin- 
orbit coupling, and superconducting proximity effects (usu- 
ally with conventional s-wave superconductors) again provide 
the key to overcoming these challenges. Note that our dis- 
cussion from Sec. |III A] on the proximity effect induced in a 
ID system applies to the 2D case with only trivial modifica- 
tions. We again stress that while the influence of a nearby 
superconductor can be crudely modeled by simply adding a 
pairing term to the Hamiltonian for the 2D electron system, 
at this level one misses important physics. Hybridization with 
the parent superconductor also effectively redu ces the energy 
scales intrinsic to the 2D system^^ see Sec 



B. 



5/2 fractional quantum Hall effect 



III A 



pendix|A| Both effects are important to keep in mind for all of 
the engineered heterostructures discussed below. 



Although our focus is on reviewing new routes to Majorana 
fermions, we would be remiss to not at least briefly discuss 
the fractional quantum Hall state observed^^ in GaAs quan- 
tum wells at filling factor u = 5/ 2 — which remains a leading 
candidate in this search. The basic question we would like to 
explore is how a 2D electron gas (2DEG) subjected to a strong 
perpendicular magnetic field can realize the physics of a topo- 
logical 2D spinless p -\- ip superconductor. The connection 
between these very different systems was first elucidated in 
highly influential work by Read and Green.'^ 

Because the field quenches the kinetic energy and induces a 
Zeeman splitting, obtaining a spinless regime in such a 2DEG 
is rather natural (though not guaranteed). The onset of p -\- ip 
'superconductivity', by contrast, is far subtler. To see how this 
arises we will first examine a 2DEG with a half-filled lowest 
Landau level (u = 1/2). Assuming perfect spin polarization, 
the system may be described by the Euclidean action 



d rdrip 



dr 



(-iV - f A)^ 
2m 



(75) 



where V is a 'spinless' fermion operator for electrons with 
mass m, V X A = Bz is the applied field, and ^int encodes 
Coulomb interactions. The physics of = 1/2 can be el- 
egantly captured by decomposing the electron ^ in terms of 
a 'composite fermion' / bound to two $o = ^ flux quanta 
using Chem-Simons theory!^ In this framework the action 
becomes 



{dr - ia°) + 



(-iV 



2m* 



4^0 



f 
(76) 



Here m* is an effective mass and is the Chem-Simons 
field. The temporal component serves as a Lagrange multi- 
plier that pins the statistical Chern-Simons flux to the compos- 
ite fermion density via V x a = —2^opfz. Since the mean 
density of a half-filled Landau level is {f^f) = B/{2^o), on 
average the attached flux exactly cancels the applied magnetic 
field. Composite fermions then behave similarly to electrons 
at B = 0: they form a 'composite Fermi sea' and can prop- 
agate in straight lines over long distances despite the strong 
magnetic field. 

This composite Fermi sea can in principle undergo a BCS 
instability, just as for a conventional metal. Moore and Read^ 
explored the possibility of p -\- ip composite fermion pairing 
and proposed the following lowest Landau level wavefunction 
for such a state. 



MR 



= Pf 



(77) 



andAp- with Ib the magnetic length and Zj the complex coordinate 
for particle j. The {zi — zj)'^ factors (roughly) correspond to 
the attached Chem-Simons flux, whereas the Pfaffian is the 
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real-space wavefunction for p + zp-paired composite fermions 
in the weak-pairing phase^^. Universal topological properties, 
such as the existence of chiral Majorana edge states and Majo- 
rana zero-modes bound to vortices, are shared by the Moore- 
Read state and the 2D spinless p + zp superconductor explored 
earlierPS 

At u = 1/2, such a pairing instability does not arise 
experimentally — the composite Fermi sea is stable and un- 
derlies the formation of an interesting compressible 'com- 
posite Fermi liquid' phase. A compelling body of theoret- 
ical evidence^^, however, indicates that the Moore-Read state 
(or its particle-hole conjugat d^^^*^^^* ) provides an energetically 
very competitive candidate for the measured plateau in the 
half-filled second Landau level. Very likely, either the Moore- 
Read state or its particle-hole conjugate emerge as the ground 
state over some range of density, quantum well width, mobil- 
ity, etc., and a growing set of experiment^^'Ell indeed support 
this possibility. For more details on this interesting subject 
we refer readers to Read and Green and the comprehensive 
review by Nayak et al^ 



C. 'Intrinsic' p -\- ip superconductivity: Sr2Ru04 

In rare cases, p -\- ip superconductivity can emerge 'in- 
trinsically' through interactions in a material. At present 
Sr2Ru04 — a layered compound with a somewhat com- 
plex, spin-de generate Fermi surface deriving from Ru d- 
orbitalJ^^SSIl — constitutes the best experimental candidate for 
such a superconductor. While the precise nature of the su- 
perconducting state that appears below Tc = 1.5K remains 
unsettled (see, e.g., Ref. 1 1701) , a variety of experiments sup- 
port the onset of spin-triplet Cooper pairing an d spontan eous 
time-reversal symmetry breaking in this system.'^^E^^E^^Ulil 



Recall from Sec. JIB that spinful 2Dp-\-ip superconductors 



allow for 1^ half quantum vortices that bind stable Majorana 
zero-modes. In this context, the recent experiments of Jang 
et al^ are particularly fascinating. These authors employed 
torque magnetometry to measure the magnetization of annu- 
lar, mesoscopic Sr2Ru04 samples as a function of an applied 
magnetic field B. With the field oriented perpendicular to the 
layers, increasing B produced discrete jumps in the magneti- 
zation at certain field values associated with nucleation of an 
ordinary |^ vortex in the sample. Remarkably, repeating the 
same experiment in the presence of a fixed in-plane field com- 
ponent 'fractionalized' these magnetization jumps into steps 
half as large — consistent with the entry of half quantum vor- 
tices. Precisely why the in-plane field should stabilize these 
defects is presently unclear, though Ref. 79 discusses one pos- 
sible scenario. (Note that Ref. 75 proposed applying perpen- 
dicular fields to stabilize half quantum vortices.) 

A few cautionary remarks are in order regarding 
Sr2Ru04 — and likely any 'intrinsic' p-\-ip superconductor — 
as a setting for Majorana physics. First, since time-reversal 
symmetry is broken spontaneously p-\-ip and p — ip pairings 
are degenerate, and domains featuring both chiralities will 
generally exist in a given crystal (see, e.g., Refs. ll75Hl77l) . 
These domains will complicate the edge-state structure rela- 



tive to the toy model discussed in Sec. |IIB[ Second, half quan- 
tum vortices need not trap Majorana zero-modes in Sr2Ru04 
crystals consisting of N > 1 layers. Consider, for instance, 
a half quantum vortex threading a Sr2Ru04 bilayer at T = 
where phase fluctuations can be neglected. In the artificial 
limit where the layers decouple, the vortex binds one Majo- 
rana zero-mode in each layer; restoring the interlayer cou- 
pling hybridizes these modes and produces an ordinary, finite- 
energy state. For larger N a chain of Majorana modes will hy- 
bridize and broaden into a gapless 'band' in the ^ oo limit. 
Strictly speaking, for any odd N a single Majorana zero-mode 
must survive the interlayer coupling but in practice may prove 
difficult to disentangle from other low-energy modes. Even 
in a single-layer sample Majorana zero-modes are protected 
only by a 'mini-gap' in the spectrum of vortex bound states 
[Eq. ([39])], which for Sr2Ru04 falls in the milliKelvin range 
since the Fermi energy exceeds the pairing gap by orders of 
magnitude. 

As an aside, we briefly mention a clever idea proposed in 
Ref. 11781 for realizing Kitaev's ID toy model along an ordi- 
nary 1^ vortex line threading a layered spinful p-\-ip supercon- 
ductor such as Sr2Ru04. Neglecting spin-orbit interactions 
and interlayer coupling, the vortex binds a pair of Majorana 
zero-modes in each layer. One can view each pair as com- 
prising a single site in Kitaev's ID toy model (recall Fig. [2]). 
When coupling between nearby Majorana zero-modes is re- 
stored, Ref. 178 predicts that the topological phase of Kitaev's 
model emerges upon driving a supercurrent perpendicular to 
the layers. The small mini-gap associated with the vortex, 
however, still poses a challenge for such a setup. 



D. 3D topological insulators 



In Sec. HID we described how one can engineer a ID 
topological superconductor using 3D topological insulator 
nano wires. Here we turn to Fu and Kane's groundbreaking 
proposal for stabilizing 2D 'spinless' p-\-ip superconductivity 
using the surface of a macroscopic 3D topological insulator^ 
We will continue to focus on materials such as Bi2SeJ2IE32ll 
whose boundary hosts a single Dirac cone described by Eq. 
(71 ). For a surface located in the (x, y) plane, the Hamilto- 



nian reads 



3DTI 



d'^rtlj'^[-iv{d^ay - dycj'') - /i]?^. (78) 



Equation (78 ) yields band energies e± (k) = ±v|k| — /i which 
correspond to the upper and lower branches of the massless 
Dirac cone sketched in Fig.jSj^a). This band structure is ideal 
for forming a 2D topological superconducting phase. First, 
accessing a 'spinless' regime is trivial here: for any ji that re- 
sides within the material's bulk band gap there exists only a 
single Fermi surface as desired (rather than two as ordinarily 
arises due to spin degeneracy). Furthermore, since the elec- 
trons along this Fermi surface are not spin-polarized, p -\- ip 
pairing can be effectively induced using the proximity effect 
with a conventional 5-wave superconductor. 
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FIG. 8. (a) Single Dirac cone describing the surface states of a 3D 
topological insulator. Because the electron spins wind by 2ti upon 
encircling the Dirac cone, the proximity effect with a conventional 
s-wave superconductor drives the surface into a time-reversal in- 
variant relative of a topological 2D spinless p -\- ip superconduc- 
tor. Chiral Majorana edge states form at the boundary between 
superconductivity- and magnetically-gapped regions of the surface 
as shown, for example, in (b). Introducing an ^ vortex through 
a topologically superconducting portion of the surface binds a sta- 
ble Majorana zero-mode. This can be achieved either by applying a 
magnetic field, or by adjusting the phases on superconducting islands 
(hexagons) as shown in (c). 



One can see this explicitly by examining an effective 
Hamiltonian for the surface with proximity-induced spin- 
singlet pairing: 



H = H' 



3DTI 



Ha 



Ha 



dPv/\{^^^^^H.c.), 



(79) 
(80) 



which describes a superconductor with a fully gapped quasi- 
particle spectrum given by 



E±(k) = Ve±(k)2 + A2. 



(81) 



To understand the nature of this state it is instructive to per- 
form a unitary transformation that diagonalizes the kinetic en- 
ergy in H. In terms of operators V^±(k) that add electrons to 
the upper and lower half of the Dirac cone, the Hamiltonian 
can be written as 



~2 



^/;,(k)^/;,(-k)+i/.c. 



} 



(82) 



It is clear in this basis that the proximate 5-wave superconduc- 
tor mediates p -\- ip pairing for electrons at the Fermi level.^^ 
Figure [8ja) illustrates the physical origin of this effect: A can 
pair resonant electrons with momenta k and — k since they 
carry opposite spins, while the nontrivial Cooper pair angular 
momentum arises because the spins rotate by 27r upon encir- 
cling the Dirac cone. 



The relation between Eq. ( |82| ) and the toy model dis 
cussed in Sec 



II B I becomes manifest when /i resides far 
from the Dirac point. For a heavily electron-doped surface, 
for instance, one can safely project out the Dirac cone's 
lower half by sending ip- 



precisely onto the Hamiltonian for a 2D 'spinless' p ^ ip 
superconductor — albeit with a non-standard kinetic energy 
and pairing potential — in the topological weak pairing phase. 
By continuity, the surface forms a topological superconduc- 
tor for any chemical potential that does not intersect the bulk 
bands since the quasiparticle spectrum of Eq. ( [ST] ) is always 
fully gapped. In addition to being 'easy' to access in this 
sense, we emphasize this phase can also in principle per- 
sist up to relatively high temperatures. For one, the surface 
ca n pote ntially inherit the parent superconductor's full pairing 
g^p(89|94j ^Qj. ^i^g same reasons discussed in the 2D topologi- 
cal insulator context; see Sec. |IIIB[ Moreover, the topolog- 



ical phase captured here preserves time-reversal symmetr}^ 
[which is obvious in the original ijj^^^ basis but somewhat hid- 
den in Eq. ([82])]. This feature guarantees that the topological 
superco nducto r's gap enjoys immunity against non-magnetic 
disorder.l^^nnil Time-reversal symmetry breaking of some kind, 
however, is required to uncover the seeds of Majorana physics 
encoded in this state. 

Similar to the toy model of Sec. |IIB| chiral Majorana 
edge states form at the boundary between topologically 
superconducting and magnetically gapped regions of the 
surface.^ ^^^tf^^" Figure [sj^b) illustrates one possible architec- 
ture supporting such an interface. There, an 5-wave supercon- 
ductor generates topological superconductivity, while a sur- 
rounding ferromagnetic insulator imparts the surface beneath 
it with a Zeeman field that we assume cants the spins out 
of the (x, y) plane. The surface state Hamiltonian govern- 
ing the latter region then becomes H = H^bti + Hz, with 
Hz = —hJcPrip^a^ip. The Zeeman energy modifies the 
spectrum to e± (k) = ±y/ ('^|k|)^ + h^ — ji, so that the surface 
forms a magnetically gapped state when the Fermi level lies 
within the resulting field-induced gap in the Dirac cone. Inter- 
estingly, since the topological superconductor induced at the 
center retains time-reversal symmetry, the edge-state chirality 
depen ds on whe ther the ferromagnet cants the spins along +z 
or —z JISMEHxhis will be important to keep in mind when 



we discuss interferometry in Sec. V C 



Another useful way to lift time-reversal symmetry is to in- 



troduce an 1^ vortex in a topologically superconducting re- 
gion of the surface; exactly as for an ordinary spinless p ^ ip 
superconductor, a single zero-energy Majorana mode local- 
izes at the vortex core.^^ While vortices can always be in- 
duced by applying a magnetic field, Fu and Kane invented an 
alternative, more versatile method for creating and manipulat- 
ing Majorana zero-modes.^^ Figure [sjc) illustrates their pro- 
posed setup, consisting of an array of superconducting islands 
(hexagons) deposited on a 3D topological insulator surface. 
Here vortices appear when the superconducting phases on the 
islands wind by ±27r around a trijunction where three islands 
meet. For example, the pattern of phases in Fig. [SJc) traps 
two Majorana modes 71,2. Manipulating the superconduct- 
ing phases (by, say, driving currents across the islands) allows 
one to controllably create, transport, and remove vortices — all 
crucial ingredients for the topological quantum information 



processing schemes that we highlight in Sec. VI 



0. Equation (82) then maps 



One commonly noted obstacle to realizing Fu and Kane's 
proposal experimentally is that most 3D topological insu- 




lators studied to date do not actually insulate in the bulk. 
Rather, they contain a substantial concentration of gapless 
bulk carriers — which must be removed entirely for the physics 
described above to survive. (Section |V A| discusses the 
fate of Majorana modes coupled to gapless degrees of free- 
dom.) While important progress towards rectifying this issue 
has recently been reported — see, e.g., Refs. '183" and "184 — 
these carriers may ultimately prove to be a feature rather 
than a bug. A number of recent experiments observe su- 
perconductivity in metallic 3D top ologica l insulator materi- 
als upon doping or under pressure. ^^3^111 xhe nature of the 
resulting superconducting state is itself a fascinating prob- 
lem, and it has been suggested that these systems may con- 
stitute the first rea lization of a class of exotic 3D topological 
superconductors^^ 63 |8i | 82 | i 93^ 198] Hqyq we will concentrate on 
the simplest possibility wherein conventional 5-wave super- 
conductivity emerges in the bulk. Naively, it is tempting to 
conclude that employing such materials kills two birds with 
one stone: the problematic bulk carriers are gapped by Cooper 
pairing and simultaneously play the role of the proximate 5- 
wave superconductor in Fu and Kane's proposal. Upon closer 
inspection, however, the validity of this physical picture is sus- 
pect. In the metallic phase well-defined surface states need 
not exist at the Fermi energy once the chemical potential in- 
tersects the bulk bands. (Generally, the surface state pene- 
tration depth diverges due to hybridization with resonant bulk 
extended states. Even when the bulk and surface-state Fermi 
surfaces are well- separated in momentum space disorder can 
still induce hybridization.) So when the bulk becomes super- 
conducting, do Majorana modes still localize at the surface 
when a vortex is present? 



An important work by Hosur et a/."^™" showed that they can 
and provided the following appealing picture for the physics. 
Consider first the limit where the bulk superconducts but the 
Fermi level resides within the bulk band gap. Well-defined 
surface states then appear and realize a topological supercon- 
ducting phase via proximity with the bulk. A vortex line 
penetrating the material binds a pair of localized Majorana 
zero-modes, one at each end, effectively realizing Kitaev's 
toy model for a ID topological superconductor. Upon rais- 
ing the Fermi level the vortex line eventually transitions into 
the trivial strong pairing phase of Kitaev's model, at a criti- 
cal chemical potential that depends on the bulk band structure 
and the vortex orientation. The crucial point is that this tran- 
sition can occur well after the Fermi level first intersects the 
bulk bands. In doped Bi2Se3, for instance, Hosur et al. pre- 
dict that a vortex line binds Majorana zero-modes up until the 
chemical potential lies roughly 0.2 — 0.3eV above the conduc- 
tion band minimum. As with any 'intrinsic' superconductor 
where the Fermi energy greatly exceeds the pairing gap, the 
small mini-gap associated with the vortex will complicate the 
identification of the zero-modes. Nevertheless, as the mate- 
rial science of 3D topological insulators is perfected the su- 
perconducting variety of these compounds provides one very 
promising venue for the exploration of Majorana physics. 



(a)^ 



s-wave SC 



2DEG with Rashba 



FM insulator 



1 

i 



(b) 



(c) FI\A insulator 



IIIIIIJI 

s-wave SC i 





FIG. 9. (a) A 2DEG with Rashba spin-orbit coupling can effectively 
realize a topological 2D 'spinless' p -\- ip superconductor when in 
contact with a ferromagnetic insulator and conventional s-wave su- 
perconductor, (b) Band structure for the 2DEG when time-reversal 
symmetry is present (gray) and with a non-zero Zeeman field (blue) 
that opens up a 'spinless' regime. (The break in the spectrum ap- 
pears for clarity.) Alternative devices that support topological phases 
appear in (c)-(e). 



E. Conventional 2D electron systems 

Following Fu and Kane's work a number of authors pur- 
sued alternative routes to engineering 2D spinless p -\- ip su- 
perconductivity using more conventional materials. We will 
in itial ly explore the semiconductor-based proposal of Sau et 
al^^, though we note that similar results in related contexts 
appear in some earlier works^^^ Consider the architecture 
of Fig. [9ja) consisting of a Rashba-coupled, electron-doped 
semiconductor 2DEG (such as InAs or InSb) sandwiched by 
an 5-wave superconductor and ferromagnetic insulator. We 
take the 2DEG to lie in the (x, y) plane and assume that the 
ferromagnetic insulator's magnetization aligns along the z di- 
rection. The following effective Hamiltonian crudely captures 
the dynamics of electrons in the semiconductor, 



H = i^2DEG + Hz + H/\ 
V2 



H2DEG = j d\iP^ 



2m 

Hz = -hJ d^rV^V^V^ 
H^= [ d^vAi^t^^ + H.c.) 



(83) 
^4) 

(85) 
(86) 



where V^J.^. adds an electron with effective mass m and spin a 
to the 2DEG. Let us first understand the Hamiltonian i^2DEG 
describing the 2DEG's intrinsic couplings. Here a denotes 
the strength of Rashba spin-orbit interactions which favor 
orienting the electron spins within the 2DEG plane, perpen- 
dicular to their momentum. The gray curves in Fig. [9jb) il- 
lustrate the band structure obtained from i^2DEG- Note that 
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the energies resemble a Dirac cone at small k, similar to a 3D 
topological insulator surface, but whose lower half eventually 
bends upward at large momenta. Because of this property two 
Fermi surfaces appear for any /i above the bottom of the con- 
duction band — one of which we would like to eliminate. 

As we have seen several times before time-reversal sym- 
metry breaking, now through the Zeeman term Hz inherited 
from the ferromagnetic insulator, again overcomes this prob- 
lem. (The Zeeman field h arises primarily from electron tun- 
neling between the two subsystems and not the magnetic field 
emanating from the ferromagnet;^^ this justifies the neglect 
of orbital effects in the Hamiltonian above. Incidentally, for 
the ID wire proposals discussed in Sec. |IIIC| orbital effects 
from an applied magnetic field are generically subdominant to 
the Zeeman energy. Ultimately this key advantage allows for 
simpler setups employing wires compared to 2DEGs.) With 

> the 'Dirac-like' portion of the 2DEG's band structure 
near k = acquires a gap as shown in the blue curves of 
Fig. [9jb). When the Fermi level resides within this gap and 
one turns on a weak proximity-induced A from Eq. ([86]), the 
resulting 'spinless' metal enters a topological superconduct- 
ing phase supporting chiral Majorana edge states and Majo- 
rana zero-modes localized at vortex cores. The physics here 
is nearly identical to the 3D topological insulator proposal 
reviewed previously: the 5- wave pair field mediates p -\- ip 
Cooper pairing because Rashba coupling causes the in-plane 
spin co mpone nts to wind by 27r upon encircling the Fermi 
surface.^^^^^^(Reference l2Q5l provides an explicit mapping to 
a spinless 2D p ^ ip superconductor.) A quantitative analysis 
analogous to that carried out for ID wires i n Sec. 
that the topological phase appears providec^^ 



III C reveals 



h > \/A^ + /i^ (topological criterion), 



(87) 



which is the same criterion given in Eq. ( [68) ). The phase dia- 
gram is thus again given by Fig.[6|;c). 

Our discussion regarding optimization of the ID wire pro- 



posal from Sec. |IIIC [ applies to the present case with almost no 
modification, so here we will simply highlight a few important 
points. For concreteness suppose that the parent superconduc- 
tor in the device of Fig. ^sl) exhibits a pairing gap of a 
few Kelvin while the ferromagnetic insulator's Curie tempera- 
ture is of order lOOK. With these rough energy scales in mind 
it is interesting to ask how large the gap for the topological 
superconductor formed in the 2DEG can be in principle. The 
answer depends on the spin-orbit energy Eso = \ma^ and 
varies from a small fraction of A^c when Ego/^sc ^ 1 up 
to the full value of Age in the opposite limit Ego/ Age ^ 1, 
highlighting the importance of sizable spin-orbit coupling .^Sl 
(Capturing this physics requires a more rigorous treatment of 



the proximity effect as discussed in Sec. Ill A and Appendix 
[A[ ) Large values of Ego also endow the topological phase with 
some robustness against disorder despite Anderson's theorem 
not applying here. In the limit h/Ego <C 1 electrons at the 
Fermi surface in the 'spinless' regime are weakly perturbed by 
the induced Zeeman field, thus suppressing the pair-breaking 
effect of disorder. Operating in this limit, however, may be 
neither possible nor desirable due to competing physics. Sat- 



isfying the topological criterion in Eq. (87) requires that h 



exceed the inherited pairing field A, and in practice it will 
likely prove advantageous to engineer the interface with the 
ferromagnetic insulator such that h A. Though this might 
reduce the topological phase's gap below its theoretical max- 
imum, large Zeeman fields facilitate tuning of the chemical 
potential into the 'spinless' regime while simultaneously in- 
creasing the tolerance to long-range potential fluctuations. 

Despite the exceptional purity with which they can be 
fabricated GaAs quantum wells are (unfortunately) unsuit- 
able 2DEG candidates. Their Rashba spin-orbit energy Ego 
falls in the milliKelvin range^^^^^^ due to the lightness of 
the constituent elements; consequently disorder is almost 
certain to dominate at the extraordinarily low densities re- 
quired to access a topological phase in this material. More 
promising are In As quantum wells, which feature Ego val- 
ues of a few tenths of a Kelvin^^^ and also contact well with 
superconductor^^. The spin-orbit energy can be increased 
further still by employing heavier materials such as InSb, or 
hole-doped semiconductors which generally exhibit larger ef- 
fective masses and spin-orbit coupling strengths compared to 
their electron-doped analogs. Interestingly, in the latter sys- 
tems the 'heavy holes' can be driven into a topological / + if- 
paired state that also supports Majorana modes. (The ex- 
tra Cooper pair angular momentum arises because as one tra- 
verses the Fermi surface in a heavy hole band the spins wind 
by Gtt rather than 27r.) 

Apart from materials considerations, experimentally 
demonstrating a topological phase in the proposed structure 
in Fig.[9|;a) poses several other challenges. Forming two high 
quality interfaces — one on each side of the 2DEG — presents 
a nontrivial fabrication problem. Furthermore, the device ex- 
hibits limited tunability: whereas the Zeeman field in the ID 
wire proposal from Sec. IIIC can be easily varied, h is now 
largely fixed once the structure is fabricated. And finally, the 
electrons in the semiconductor are effectively buried in the 
sandwich structure, making it difficult to manipulate or probe 
these carriers. It may, fortunately, be possible to alleviate 
some of these challenges by employing modified setups. 

The interdigitated ferromagnet- superconductor device 
shown in Fig. ^c), for example, can realize a topological 
phase while requiring lithography on only one side of the 
2DEG.^^^ In such a structure the electrons inherit periodically 
modulated Zeeman and pairing fields due to the interdigi- 
tation. Provided their Fermi wavelength exceeds the finger 
spacing, however, the electrons effectively feel the spatial 
average of these quantities, resulting in a robust topological 
phase under similar conditions to the sandwich structure in 
Fig.[9|;a). (A topological phase can appear even in the oppo- 
site limit, though with a diminished gap.) Somewhat more 
elaborate interdigitated setups also allow one to electrically 
generate vortices binding Majorana zero-modes by applying 
currents to modulate the superconducting phases along the 
fingers. 

It is even possible to eliminate the ferromagnetic insulator 
altogether and drive a transition into the topological phase us- 
ing an external magnetic field, similar to ID wire setups. As 
alluded to earlier the principal reason for the ferromagnetic 
insulator in Figs.[9|;a) and (c) is that one wishes to generate a 
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Zeeman field that cants the spins out of the 2DEG plane while 
avoiding orbital effects that would accompany an applied per- 
pendicular magnetic field. In-plane fields largely circumvent 
unwanted orbital effects but unfortunately do not open a 'spin- 
less' regime — at least in a semiconductor with only Rashba 
coupling. Reference 1205 1 showed that in-plane fields can gen- 
erate topological superconductivity provided one employs a 
2DEG grown along the (110) direction with strong Rashba 
and Dresselhau^^ spin-orbit coupling (such as InSb); see 
Fig. [9jd). The special feature of this growth direction is that 
these two kinds of spin-orbit interactions conspire to rotate the 
plane in which the spins orient away from the 2DEG plane, so 
that an in-plane field plays a similar role to the ferromagnetic 
insulator in a Rashba-only setup. Another system that may 
eliminate the need for a ferromagnetic insulator is NaCo02. 
Recent first-principles calculations predict that this material is 
a (conventional) bulk insulator with surface states exhibiting 
very strong Rashba coupling and spontaneous time-reversal 
symmetry breaking that opens a broad 'spinless' regime. ^"^^ 
These predictions would be interesting to explore experimen- 
tally using ARPES. 

Early on Lee suggested another route to engineering topo- 
logical superconductivity, using a 2D fully spin-polarized fer- 
romagnetic metal adjacent to a spinful bulk p -\- ip supercon- 
ductor such as Sr2Ru04 [Fig.|9je)].^^^ Essentially, the ferro- 
magnetic metal 'filters out' one spin component from the par- 
ent superconductor and realizes a 2D spinless p-\-ip supercon- 
ductor due to the proximity effect. Later it was realized that 
even a conventional 5-wave superconductor can drive a fer- 
romagnetic metal into a topological phase provided apprecia- 
ble spin-orbit coupling appears at the interface.^ One virtue 
of these setups is that the topological phase can in principle 
exist without requiring fine-tuning of the chemical potential, 
though finding suitable 2D ferromagnets poses an experimen- 
tal challenge. A more exotic alternative proposal predicts 
that topological superconductivity can appear in proximity- 
coupled systems realizing a novel 'quantum anomalous Hall' 
state.^^^ This possibility provides strong motivation for ex- 
perimentally pursuing this as yet undiscovered phase of mat- 
ter. 



V. EXPERIMENTAL DETECTION SCHEMES 

One very simple (but also highly indirect) way of inferring 
the existence of a Majorana mode is through the detection of 
a topological quantum phase transition in the bulk of a ID or 
2D superconductor. Consider for example the ID spin-orbit- 
coupled wire proposal reviewed in Sec. [Ill C| where one can 
tune between trivial and topological superconducting states 
simply by turning on a magnetic field. Observing the bulk 
gap collapse and then reopen as the field strength increases 
would provide strong evidence for the onset of topological 
superconductivity and, by extension, the appearance of Ma- 
jorana end states in the wire. A topological phase transition 
in this setting also manifests itself in thermal and electrical 
transport 
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FIG. 10. A ID spin-orbit-coupled wire draped across an insulator-s- 
wave superconductor junction allows one to detect Majorana modes 
via transport. In (a) the wire's left half is gated into a 'spinless' metal- 
lic regime while the right half forms a trivial gapped superconduc- 
tor; in this setup the zero-bias conductance of the junction vanishes. 
When the right half of the wire instead forms a topological phase as 
in (b), hybridization between the 'spinless' metal and the Majorana 
7i produces a quantized 2e^ /h zero-bias conductance. The same 
results apply when the wire's left half forms a spinful metal, pro- 
vided repulsive interactions are present. This geometry can be read- 
ily adapted to detect isolated Majorana modes in many other settings. 



desirable. Below we review three classes of such measure- 
ment schemes, based on tunneling, Josephson effects, and in- 
terferometry. 



A. Tunneling signatures of Majorana modes 

Tunneling spectroscopy provides a powerful and concep- 
tually appealing method for detecting Majorana zero-modes. 
To illustrate the physics we will focus on the experimentally 
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as well as Coulomb blockade experiments 
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Clearly, however, more direct probes of Majorana modes are 



accessible geometry shown in Figs. [TO[a) and (b), where a 
long spin-orbit-coupled wire subjected to a magnetic field ex- 
tends across an insulator-5-wave superconductor junction 
In both cases we assume that the wire's left half is gated into 
a 'spinless' regime and remains metallic. The proximity ef- 
fect with the 5-wave superconductor, however, drives the right 
half into a trivial gapped state in (a) but a topological phase 
supporting Majorana modes 71,2 in (b). (One can tune be- 
tween these configurations by adjusting the chemical poten- 
tial in the right half.) We would like to contrast the conduc- 
tance of these two setups — particularly at zero bias — when 
current flows from the 'spinless' metal into the superconduc- 
tor. Our approach will follow closely Refs. .220i and especially 
122 II which emphasizes universal features of the problem. 

We initially attack the trivial case of Fig. [TOja). The only 
available low-energy degrees of freedom reside in the metal- 
lic region, which is taken to lie along x < 0. We there- 
fore model the system by an effective Hamiltonian H = 
Hmetai + ^junction in which the first term describes the lin- 
earized kinetic energy near the Fermi level at x < while 
the second captures terms at x = imposed by the supercon- 
ductor. Defining spinless fermion operators ipR/L describing 
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right/left-moving excitations in the wire, we have 



metal 



(88) 



with vf the Fermi velocity. One can rewrite i^metai in terms 
of a single chiral field ^|J{x) defined over all x such that > 

0) = i^L{—x) and iIj{x < 0) = ^Ijr{x), yielding 



metal 



dx (^—ivFt/j^dxip) 



(89) 



The leading local potential and pairing terms induced at the 
junction read 



dx [mp^ip + A{iljdxilj + H.c.)] S{x) 



(90) 

In particular, the A term (which requires a derivative by Fermi 
statistics) allows Cooper pairs to hop from the metallic region 
of the wire into the superconductor. 

The full Hamiltonian H takes the form of a Bogoliubov-de 
Gennes equation and is therefore diagonalized with quasipar- 
ticle operators carrying energy E of the form 



-f 

J — c 



dxe 



-[Pe{x)iP{x)^He{x)iP^{x)], (91) 



where Pe and He determine the particle- and hole-like am- 
plitudes for the quasiparticle wavef unctions. The conductance 
we are after follows from the scattering matrix S{E) that re- 
lates states incident on the superconductor to the reflected 
states. Recalling the definition of in terms of left and 
right movers, the S'-matrix components are defined by 



Spp{E) Sph{E) 



Pe{-oc) 
He{-^) 



(.92) 



Since \Sph{E)\'^ is the probability that an incident electron 
at energy E Andreev reflects as a hole at the junction, passing 
charge 2e into the superconductor, a bias voltage V applied 
across the junction generates a current 



_ ^ r 



dE\SpH{E)\' 



The differential conductance G{V) = ^ follows as 

G{V) = ^\SpH{eV)\'. 



(93) 



(94) 



Very general arguments adapted from Refs. "179" and 180 
strongly constrain the iS-matrix at zero energy, and hence the 
zero-bias conductance. First, particle-hole symmetry of the 
Bogoliubov-de Gennes equation dictates that Te = — r^, 
which in turn implies S{E) = cr^ S*{—E)cr^ . This relation, 
together with unitarity of the ^-matrix as required by current 
conservation, restricts 5'(0) to one of two possible forms: 



5(0) 



e*" 
e-^" 



e'^ 




(95) 



for some phases a^^. The purely diagonal case corresponds 
to the onset of perfect normal reflection — with unit probabil- 
ity an electron reflects as an electron and similarly for holes — 
and hence a vanishing zero-bias conductance. In contrast, the 
off-diagonal possibility yields perfect Andreev reflection; here 
electrons scatter perfectly into holes and vice versa, yielding 
a quantized zero-bias conductance G = 2e^ /h. These very 
different limits represent renormalization group fixed points 
at which the superconductor imposes either perfect normal re- 
flecting or perfect Andreev reflecting boundary conditions on 
the metal at low energies. 1^ We stress that Eq. ( 95 ) relies only 
on the metal being 'spinless', and holds even when Majorana 
zero-modes are present as in Fig.p^b). 



By explicitly calculating the S'-matrix for Eqs. (89) and 



(90) it is straightforward to show that 5(0) is purely diago- 
nal so that the zero-bias conductance vanishes in Fig. [TOja). 
The A term in injunction does permit Cooper pairs to tunnel 
into the superconductor, but because the metal is 'spinless' 
the probability vanishes at zero energy due to Pauli block- 
ing. In other words the non-interacting 'spinless' metal we 
have treated so far generically flows at low energies to the per- 
fect normal reflection fixed point. It is, however, essential to 
understand the impact of interactions — which transform the 
metal into a Luttinger liquid — on this result. Reference l221l 
demonstrates via bosonization that this fixed point remains 
stable even in the interacting case except when the wire ex- 
hibits very strong attractive interactions. Remarkably, for such 
an attractive wire a pair of asymptotically decoupled Majo- 
rana modes emerges dynamically at the junction and drives 
the system to the perfect Andreev reflection fixed point with 
quantized conductance.!^ 

Dramatically different physics arises for the topological 



setup in Fig. [TO[b) due to the Majorana modes. We will as- 
sume zero overlap between 72 and 71 and now focus on a 
junction Hamiltonian 



junction 



J — ( 



<^^7i('0 — iIj)5{x) 



(96) 



that hybridizes 71 and the 'spinless' metal. [The u and A 
terms from Eq. ( [9Q| ) are qualitatively unimportant so we ne- 
glect them for simplicity.] For any t 7^ the Majorana 71 no 
longer represents a zero-energy mode since [injunction, 7i] 7^ 
0. A single Majorana zero-mode can never exist on its own, 
however, so where is 72 's partner? The answer is that the 
zero-mode described by 71 when t = gets absorbed into 
the metal where it becomes a delocalized plane wave; ex- 
plicitly, one can readily verify that for t ^ 72's partner is 
7i dx{\lj + ijj^)}'^'^ The hybridization with 71 lead- 

ing to this delocalized Majorana plane- wave mode mediates 
perfect Andreev reflection at the junction at low energies. In- 
deed, extracting the ^-matrix from the Hamiltonian given by 



Eqs. ( |89) ) and ( [96| ) yields a Lorentzian conductance 

2e2 
h 



GiV) 



1 



eVvp 
2^2 



(97) 



that collapses to 2e^ /h in the zero-bias limit. From a renor- 
malization group perspective, coupling to 71 drives the non- 
interacting 'spinless' metal to the perfect Andreev reflection 
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fixed point characterized by a quantized zero-bias conduc- 
tance, in stark contrast to the trivial setup of Fig.[TOja). This 
result survives even for a Luttinger liquid unless strongly re- 
pulsive interactions with a Luttinger parameter g < 1 /2 are 
present; there the repulsion obliterates the zero-bias conduc- 
tance peak entirelyP2II 

What happens when the metallic region of the wires in 
Fig. [To] carry spin? In this case the arguments leading to Eq. 
( [95] ), which underlies the conductance dichotomy for the se- 
tups in (a) and (b), provide much weaker constraints on the 
S'-matrix. Furthermore, a local pairing term at the junction 
Ao[ijj^{x = 0)V^|(x = 0) + H.C.] now allows singlet Cooper 
pairs to tunnel into the superconductor without Pauli blocking 
at zero energy. When the spinful metal impinges on a topolog- 

;b) a quantized 2e^ /h zero- 
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ical superconductor as in Fig. 
bias conductance nevertheless emerges from coupling to the 
Majorana 71 , both at the non-interacting level and over a range 
of interactions. For a non-interacting spinful metal adjacent 
to a trivial superconductor as in Fig.[TOja), however, Aq pro- 
duces an iS-matrix that yields a zero-bias conductance rang- 
ing from to Ae^ /h depending on parameters — potentially 
making it difficult to differentiate with the topological case. 
Fortunately this result is non-generic. Arbitrarily weak repul- 
sive interactions drive the system to the perfect normal reflec- 
tion fixed point at low energies, leading to a vanishing zero- 
bias conductance just as for a 'spinless' metal. Thus sharp 
tunneling signatures of Majorana zero-modes appear also in 
the spinful case. Adding more channels, however, obscures 
these signatures since the conductance in both the trivial and 
topological setups is then non-universal, at least in the free- 
fermion limit.'^ 

While our discussion so far centered on a specific geome- 
try involving ID wires, the conclusions apply far more gen- 
erally. Nowhere in this analysis did we use the fact that the 
Majorana 71 derived from a topological region of a ID wire. 
In fact the quantized 'zero-bias anomaly' captured above has 
been discuss ed in numer ous contexts from several different 
perspectives ^^SElEIIEllEISI^j^^ ^ general property of spin- 
less or spinful ID metals tunneling onto an isolated Majo- 
rana mode. The setups of Fig. 10 can be readily adapted 
to probe Majorana modes in 2D topological insulator edges, 
half-quantum vortices in a 2D spinful p -\- ip superconductor, 
chiral Majorana edge states^^^, etc. (Though one should keep 
in mind that the small mini-gap typically associated with vor- 
tices and edge states can place stringent limits on temperature 
and resolution.) 

One might view the tunneling signatures of Majorana 
modes discussed above as somewhat less than a 'smoking 
gun' detection method since zero-bias anomalies can arise 
from unrelated sources in mesoscopic systems. We are some- 
what sympathetic to this perspective but stress the following 
points. First, whereas Majorana modes produce a zero-bias 
conductance of 2e^ /h this quantized value need not appear 
when tunneling into a conventional 'accidental' low-energy 
mode. A second, more important point is that in many of the 
proposals we reviewed this conductance peak can be control- 
lably brought in and out of resonance by inducing a transi- 
tion between topological and trivial phases (by, say, adjusting 



the magnetic field or gating). While measuring a zero-bias 
peak may not by itself constitute a smoking-gun signature of 
a Majorana mode, observing this peak collapse and revive in 
accordance with theoretical expectations arguably would. 

Majorana zero-modes provide several other note- 
wort hy transpo rt signatures, including through current 
noise. EniI13[226] particular, coupling a pair of ID metallic 
wires to a topological superconductor introduces Majorana- 
mediated 'crossed Andreev reflection'; such processes arise 
when a Cooper pair enters the superconductor by combining 
a single electron from each metal and generate maximally 
correlated current noise in the two wires. ^^^^^^ The tun- 
neling conductance through a quantum dot has also been 
predicted to change qualitatively when the dot couples to 
a Majorana mode,^^^^^^ and a related setup allows one to 
probe the lifetime of the ordinary fermionic state formed by 
a pair of Majorana modes in the presence of 'quasiparticle 
poisoning Additional transport signatures arise when a 
Majorana mode couples to a nanomechanical resonator.ESl 
Spin-polarized scanning tunneling microscopy may also be 
used to identify fingerprints of Majoranas.^^^ Finally, in 
a mesoscopic topological superconductor where charging 
energy is important, spatially separated Majorana modes can 
mediate non-l ocal ele ctron 'teleportation' that can be detected 
with transport.EUSa 



B. Fractional Josephson effects 

Recall from the introduction that a system with 2A^ well- 
separated Majorana zero-modes 71,..., 2 at exhibits 2^ degener- 
ate ground states. By defining conventional fermion and num- 
ber operators 



/. = (7: 



•n2j)/2, 



(98) 



the ground-state manifold can be conveniently labeled by 
|ni,...,nAr), where Uj = 0,1 specify topologically pro- 



tected qubit states. [The pairing of Majoranas in Eq. (98) 
is completely arbitrary but always sufficient; ground states 
labeled with different pairings are simply connected by uni- 
tary transformations.] The experimental detection methods 
discussed so far allow one to deduce the existence of Majo- 
rana modes but provide no information about the qubits they 
encode. One way to extract this information is to prepare the 
system into a ground state and then adiabatically bring two 
Majorana modes — say 71 and 72 — in close proximity so that 
their wavefunctions overlap appreciably. The resulting hy- 
bridization of these modes can be modeled by a Hamiltonian 



Sl\l2 = e(ni — 1/2). Taking e > for concreteness. 



the system remains in a ground state if ni = whereas the 
fusion of 7i and 72 yields an extra finite-energy quasiparticle 
if ni = 1. One can thus read out the state of ni by detect- 
ing the presence or absence of such a quasiparticle. (Note that 
multiple measurements may be required since the system can 
form a superposition of ni = and 1 states.) Fusing Majo- 
ranas across a Josephson junction both enables qubit readout 
along these lines"^^ and provides an unambiguous fingerprint 
of these modes. This approach is particularly well- suited for 
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FIG. 11. (a) Basic experimental setup required to observe the frac- 
tional Josephson effect stemming from Majorana modes fused across 
a superconductor-insulator-superconductor junction. Purple regions 
indicate s-wave superconductors (with phases (t)L,R) that drive the 
green regions into a ID topological state; dashed regions are trivially 
gapped. The Majoranas 71,2 mediate a component of the Josephson 
current that is 47r periodic 'm(j)R — (j)L. When the barrier in the junc- 
tion is replaced by a superconductor with phase (j)M as in (b), 71,2 
mediate a second type of unconventional current that is An periodic 
in and can be isolated with Shapiro step measurements. The 
setup in (c), while superficially similar to (a) and (b), yields only 
conventional Josephson physics with 27r periodicity in(j)L,R. 



ID topological superconductors, so we will confine our dis- 
cussion to this class of systems. 



Figure [TT] sketches the basic type of setup required. Here 
two ID topologically superconducting regions emerge due to 
the proximity effect with 5-wave superconductors that are sep- 
arated by an insulating barrier, while dashed regions are as- 
sumed to be trivially gapped. The topological segments can 
arise from any number of systems — a 2D topological insula- 
tor edge, ID spin-orbit-coupled wires, 3D topological insula- 
tor nano wires, etc. — but should be long enough that the outer 
Majoranas 73,4 overlap negligibly with the central Majoranas 
7i,2-^^^ The insulating barrier, however, should be sufficiently 
narrow that 71 and 72 hybridize strongly. Given this setup, 
our objective is to understand the zero-bias current / flowing 
across the Josephson junction as one varies the phase differ- 
ence /S.(j) = — (j)L between the right and left 5- wave super- 
conductors. 

This current consists of two contributions, / = + /e- 
The first, /2e, denotes the conventional Josephson current that 
arises from Cooper-pair tunneling across the insulating barrier 
and is 27r-periodic in A^. As originally shown by Kitaev^ the 
hybridized Majorana modes 71^2 mediate a new contribution 
le which is our primary focus. To extract the salient univer- 
sal features of /g in a very direct way we will model the two 
topological regions of Fig. [TTJa) as A^-site chains described 
by Kitaev's toy lattice model given in Eq. ([2]). Furthermore, 
we fine-tune /i = and t = A for each region so that the 
Majorana zero-modes can be trivially identified. Defining op- 
erators c^i^i^^ that add fermions to the left/right topological 
segments, the full Hamiltonian for the junction is then taken 



H= ^a + ^r, (99) 

a=L/R 
t ^"^ 

(101) 



where Hy describes single-electron tunneling across the bar- 
rier with strength F > 0. 



Let us recall the following two facts derived in Sec. II A 
(i) Ha supports Majorana zero-modes localized at sites 1 and 
A^ of each chain [see Fig.[2jb)] and (n) the zero-modes 71,2 
at the junction are related to the lattice fermion operators by 
CLN = e-^^^/2(7i+i7;)/2andcm = 6-^^^/2(7^ + ^72)72. 
Here 7^ and 72 hybridize with Majorana fermions at neighbor- 
ing sites and form conventional finite-energy fermions. One 
can therefore project H onto the zero-energy subspace of Ha 
by sending 



1 

cln 



71, cri 



-ic/)R/2 



72, 



(102) 



which yields an effective low-energy Hamiltonian 

i^eff = ~ 2 "^"^^ V ^ / ^^^^^ 

= -rcos( (ni-1/2). 



(103) 



Crucially, the occupation number ni is a conserved quantity 
since [i^eff, ^1] =0. Thus if the system begins in a state with 
rii = n\, then varying the phase difference A^ across the 
junction yields a Majorana-mediated current 



2e (i^eff) 
h dAcj) 



eV /A^ 

— r sm 

2h \ 2 



(2nl-l). 



(104) 



(If the system does not form an ni eigenstate, then the current 
le above emerges with a probability determined by the relative 
amplitude forn^ =0,1 states.) 

Equation (104) reflects Si fractional Josephson effect — /g 



originates from tunneling 'half of a Cooper pair across the 
junction and exhibits An periodicity in A0. The first prop- 
erty is easy to understand: Cooper-pair hopping dominates 
the Josephson current in conventional 5-wave superconduc- 
tors because the bulk gap suppresses single-electron tunnel- 
ing, but since Hr couples zero-modes this suppression disap- 
pears. The 47r periodicity is much subtler. Indeed, the original 
Hamiltonian in Eq. ( 99 ) is clearly 27r periodic in both (j) l and 
(})R, SO how can the current /g exhibit 47r periodicity? This 
is possible because while the Hamiltonian is 27r periodic, the 
physical states are not. Suppose, for example, that A^ = 
and n\ = 1 so t hat the system begins in a ground state of 
i^eff in Eq. (103) with energy Ei = — F/2. Because ni is 
conserved, after advancing A^ by 27r the system ends in a 
physically distinct excited state with energy Ef = +F/2 and 
hence an extra finite-energy quasiparticle at the junction .ESI 
Global fermion parity conservation dictates that the system 
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can decay back to a ground state only if fermions can transfer 
between the inner and outer Majoranas in Fig. [TTJa), which 
we precluded due to their large spatial separation. Advancing 
A(/) by 27r a second time restores the ground state with which 
we started. This doubled A ^-periodicity in the physical states 
underlies the fractional Josephson effect uncovered above. 

Measuring a 47r-periodic contribution to the Josephson cur- 
rent would undoubtedly qualify as a smoking-gun signature of 
Majorana fermions. Moreover, because the sign of the current 
le in Eq. ( 104 ) is tied to the occupation number n\ of the Ma- 
joranas at the junction, this technique enables qubit readout as 
claimed above."^^ The following experimental realities should, 
however, be kept in mind. First, /g must be disentangled from 
the (potentially much larger) 27r-periodic component /2e that 
flows in parallel. To obtain a crude order-of-magnitude esti- 
mate for le, with F ~ IK in Eq. (104) the associated crit- 
ical current is — eT/(2h) ~ lOnA, which roughly sets 
the required current resolution. Second, due to the finite ex- 
tent of their wavefunctions the outer Majorana modes 73,4 of 



Fig. 11 a) will inevitably couple to 71,2 with a characteristic 
energy 5E oc e~^^^ (L denotes the size of the topological 
regions and ^ is the topological phase's coherence length). 
Though exponentially suppressed in L/^, this hybridization 
spoils con serva tion of ni and restores 27r periodicity of the 
current I^^^^To circumvent this problem A^ should cycle 
on a time scale that is short compared to hV/SE'^P^ but long 
on the scale set by the inverse bulk gap. Third, inelastic pro- 
cesses involving stray quasiparticles — which can appear, e.g., 
because the system was imperfectly initialized or due to ther- 
mal excitation — provide another means of switching the value 
of ni to relax back to the ground state of i^eff. This, too, 
can restore 27r periodicity of /g if A0 cycles on scales much 
longer than the typical switching time; on shorter scales /g 
will exhibit 'telegraph noise' where the sign of the Majorana- 
mediated current abruptly changes with time, which would 
also be remarkable to observe.^^ Even in the presence of re- 
laxation processes signatures of the fused Majorana modes at 
the junction ap pear through the current noise spectrum^^ 
and in transientP^H, 

A second type of fractional Josephson effect arises when 
the barrier material in the junction forms a superconductor 

b).l^ The influence of 
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with phase as shown in Fig. 
the superconducting barrier can be crudely modeled by sup- 
plementing the Hamiltonian in Eq. ([99]) with a new term 

Hm = AMie'^^CLNCRi + H.c.) (105) 
that Cooper pairs fermions at the inner ends of the topologi- 



cal regions. Using Eq. (102) to project the full Hamiltonian 
onto the low-energy subspace formed by the Majoranas then 
produces a modified effective Hamiltonian 



i^eff i^eff + Am cos 



(m - 1/?|06) 



The meaning of the Am term above can be simply understood 
by promoting e*"^^'^'^ to Cooper-pair creation operators; Am 
then clearly reflects processes whereby a Cooper pair in the 
central region fractionalizes, with half entering the left topo- 
logical segment and the other half entering the right. This 



Cooper-pair splintering allows for a current 

L —-. — sm I (t)M \2n\ - 1) (107) 



injected into the middle region to be carried away in equal 
parts into the left and right topological superconductors. Like 
the fractional Josephson current in Eq. (104), exhibits 47r 
periodicity in both (j)L and (j)R because of conservation of ni, 
and also allows one to read out the initial occupation number 
n\ . A promising feature of this setup is that Shapiro step mea- 
surements can be used to isolate the I'^ contribution from the 
parallel components /2e and /g-^^^ 

Although the Majorana-mediated currents /g and were 
derived in a fine-tuned toy model, their anomalous 47r 
periodicity has a topological origin and thus arises far 
more generally (subject to the caveats noted above). In- 
deed, fractional Josephson effects have been capt ured in 
more realistic models for ID p-wave superconductor^^^, 
2D topological insulator edge s^^*^^^', 3D topological insula- 
tor surfaces49 242^ single^ ^SEmSEIl^j^^ multi-c hannel ^ 

spin-orbit-coupled wires, among other systems^^^^"^^ and 

have even been shown to survive when capacitive charging 

WE 



energy is incorporated^. We also note that a long Josepshon 
junction formed by 2D topological superconductors supports 
an unconventional Fraunhofer pattern arising from chiral Ma- 
jorana edge states. 1^ 

Some experimental setups more easily lend themselves to 
observing fractional Josephson effects than others since this 
phenomenon requires stabilizing extended topological regions 
on both sides of the junction as shown in Figs.[TTJa) and (b). 
For instance, topological insulator edges may be ideally suited 
for this type of experiment since there the topological phase 
is 'easy' to access in the sense discussed in Sec. |IIIB| We 



should caution that Josephson physics in the setup of Fig. 
[TTJc) — which in realizations such as ID spin-orbit-coupled 
wires should be much simpler to realize than those of Figs. 
[TTJa) and (b) — is always 27r periodic in the superconducting 
phases. Perhaps the simplest way to see this is by observing 
that in the limit where 71 and 72 hybridize significantly across 
the junction in Fig.[TTJc), the central region is in no meaning- 
ful way topological. 



C. Interferometry 

Several interferometric schemes have been pro- 
posed for det ecting Major ana modes in topological 
superconductors.' ^^ ' ^^ ' ^^^ ' ^ ^^'^^^Jt^ While the setups required 
generally pose greater fabrication challenges compared to 
the measurement techniques reviewed earlier, the signatures 
that appear are unambiguous and conceptually illuminating. 
For concreteness we will concentrate on interferometers 
employing 3D topological insulators; we stress, however, that 
similar ideas can be applied to many other realizations of 2D 
topological superconductivity reviewed in SecJlV| 

Consider the geometry shown in Fig. |12|^a) where an 
5-wave superconductor and two ferromagnetic insulators — 
importantly, with opposite magnetizations — reside on a 3D 
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(b) 
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FIG. 12. Interferometers fabricated from ferromagnetic insulators 
and s-wave superconductors deposited on a 3D topological insulator 
surface. In (a) the ferromagnets have opposite magnetizations, yield- 
ing a conventional chiral edge state (double arrows) that fractional- 
izes into chiral Majorana modes (single arrows) around the super- 
conductor. With a flux <^ — threading the center, the zero-bias 
conductance describing current flow from the left domain wall into 

9 2 

the superconductor vanishes for even Nv and is quantized at for 
odd Nv. In (b), when $ = Josephson vortices flowing through 
the interferometer produce a vortex current that oscillates with the 
charge Q on the central island (due to the Aharonov-Casher effect). 
When $ = 1^ , however, the oscillations disappear as a consequence 
of non-Abelian statistics. 



topological insulator surface P^^^ The ferromagnets drive 
the upper and lower portions of the surface into gapped quan- 
tum Hall states with Hall conductivities = ±^^^^^ 

2 

Since the Hall conductivities differ by ^, an 'ordinary' chiral 
edge state (denoted by double arrows in the figure) appears 



at each magnetic domain wall. As described in Sec. IV D 
the s-wave superconductor drives the surface beneath it into 
a time-reversal-invariant 2D topological superconductor sup- 
porting chiral Majorana edge states whose chirality follows 
from the neighboring ferromagnet's magnetization. When the 
ordinary chiral mode at the left magnetic domain wall meets 
the superconductor, it therefore fractionalizes into a pair of 
co-propagating Majorana niodes and then recombines at the 
right magnetic domain wall.EMini 

Of interest here is the conductance characterizing current 
flow from the left magnetic domain wall into the supercon- 
ductor, when a magnetic flux ^ pierces the central region 
of Fig. [T2ja). Following the scattering analysis of Sec. |V A 



the conductance at a bias voltage V is given by G{V) 
'^\Sph{^V)\^ , where Sph is now the amplitude for an elec- 
tron incident from the left to transmit a Cooper pair to the su- 
perconductor and exit as a hole at the right domain wall. As 
in our treatment of a spinless metal impinging on a supercon- 
ductor, particle-hole symmetry and current conservation again 
constrain the zero-bias conductance to G(0) = or In 
other words, at zero energy an incident electron first splin- 
ters into Majoranas 7^/5 localized at the top/bottom edges of 
the superconductor and then emerges with unit probability as 
either an outgoing electron c\^^ or hole Cout (superpositions 
are forbidden). 

When the flux vanishes one can deduce the conductance by 
adiabaticall y de forming the area of the superconducting re- 
gion in Fig. 12 a) to zero.^^^ In this limit an electron incident 



from the left is guaranteed to exit as an electron on the right, 
and by continuity the same must be true when the supercon- 
ducting region is finite. This can be summarized schemati- 
cally by the process 



^ 76 + 



(108) 



where 75 + i^t represents the intermediate fermionic state of 
the incident electron. The z ero-bias conductance therefore 
vanishes when ^ = Q.^^^MI Next, suppose that one threads 
flux <l> = 1^ to induce a single vortex in the superconduct- 
ing region. This introduces a branch cut indicating where the 
phase jumps by 27r, which we take to emanate from the core 



to the top edge in Fig. 12 a). Since 7t acquires a minus sign 
upon crossing the branch cut (recall Sec. |IIB| ), incident elec- 
trons now exit perfectly as holes. 
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76 + 'lit ^ 76 - 'lit 



^out 



(109) 



yielding a ^ conductance ! ^ ^^ 1 ^ By generalizing this pic- 
ture to flux ^ = one can see that the conductance os- 
cillates between quantized values G(0) = for even vortex 
number and G(0) = ^ for odd N^W^^ Observing 
these discrete conductance oscillations (which have also been 
captured in semiconductor-based systems^^^) would provide 
clear evidence for the chiral Majorana edge states underlying 
this remarkable result. Furthermore, concepts pioneered in 
the quantum Hall context"^^ 251 253 employed in related 
experiments to implement interferometric readout of the qubit 
states formed by vortex Majorana zero-modes 

References |^ and 12461 proposed interesting alternative 
probes of Majorana modes that rely on interferometry of vor- 
tices in the bulk of a 2D topological superconductor. Such 
'Abrikosov vortices' tend to behave classically, so we will dis- 
cuss an elegant follow-up proposal employing 'Josephson vor- 
tices' that more readily exhibit quantum phenomena.!^ Fig- 
ure 12 ;b) illustrates the desired setup, consisting of 5-wave 
superconductors and ferromagnets patterned on a 3D topo- 
logical insulator (though any spinless p -\- ip superconduc- 
tor realization will do here). In the center of the structure 
sits an island that hosts charge Q and a magnetic flux ^. As 
usual the innermost superconducting edge supports a Majo- 
rana zero-mode when the flux $ induces an |^ vortex. The 
inner and outer superconductors in Fig. [T2jb), however, re- 
alize a Josephson junction bridged by a thin ferromagnetic 
barrier. Consequently the chiral Majorana edge states at the 
interface hybridize across the junction and generally acquire 
a gap, but are not entirely inert. Remarkably, Grosfeld and 
Stern showed that Josephson vortices — at which the supercon- 
ducting phase difference across the junction locally winds by 
27T — trap a single Majorana zero-mode as in the case of an 
Abrikosov vortex.*^ 

When a Josephson vortex binding a zero-mode 7 flows 
rightward [which can be arranged by driving a perpendicu- 
lar supercurrent / as shown in Fig. 12^b)] interference of the 
upper and lower trajectories depends on both the charge Q 
and flux ^ on the center island.*^^*^ The Q dependence re- 
flects the Aharonov-Casher effecP^^I — an |^ flux encircling 
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charge Q acquires a phase = If ^ = ^ ' similar to the 
Aharonov-Bohm phase accumulated by a charge encircling 
a flux. With <l> = a Josephson-vortex current ly flowing 
through the interferometer therefore oscillates with Q accord- 
ing to 



= /^[l + Acos(^Q/e)], 



(110) 



where /f? and A denote the mean current and oscillation 
amplitude l^^l^ The oscillations in Eq. (110) can be detected 
by measuring the transverse voltage difference induced by the 



vortex flow as Q varies. A flux <l> 



he 

2e 



produces a Majo- 



rana zero-mode 7in inside of the interferometer and changes 
these results qualitatively. In this case taking the Josephson 
vortex in Fig. [T2jb) around the island then leads not only 
to an Aharonov-Casher phase, but also changes the sign of 
both 7 and 7in due to branch cut crossings. The latter effect 
causes the amplitude A to vanish^^^^ (which is rooted in 
non-Abelian statistics explored in the ne xt section), destroy- 
ing the vortex-current oscillations. ^^^S^ The striking depen- 
dence on Q and <l> is a dramatic manifestation of Majorana 
modes, in particular the exotic statistics they underpin, and 
would be fascinating to observe. 



VI. NON-ABELIAN STATISTICS AND QUANTUM 
COMPUTATION 




71 



(b) 



72 




FIG. 13. (a) A 2D topological superconductor with four well- 
separated vortices binding Majorana zero-modes. Under a clockwise 
braid of 71 and 72 these Majoranas not only interchange positions, 
but 71 also crosses a branch cut. Consequently the exchange sends 
71 —72 and 72 ^71. This simple picture due to Ivano\^^ leads 
to non-Abelian statistics as described in the text, (b) Performing 
meaningful exchanges of Majorana modes arising in ID topologi- 
cal superconductors (green) requires arranging wires into networks, 
(c) Corner junctions formed by 2D topological insulators similarly 
allow Majorana modes to be interchanged along the edges. Despite 
the absence of vortices the exchange statistics remains non-Abelian 
in these networks. 



The experimental realization of Majorana modes would 
pave the way to far-reaching technological innovations. On 
the most basic level, a set of Majorana-carrying vortices or 
domain walls non-locally encodes quantum information in 
the degenerate ground-state space, enabling immediate ap- 
plications for long-lived 'topological quantum memory'. In 
the longer term the prospect of manipulating that informa- 
tion in a manner that avoids decoherence would constitute 
an important breakthrough for quantum computation. This is 
made possible by the most coveted manifestation of Majorana 
fermions: non-Abelian statistics. 

Before turning to specific implementations it will be use- 
ful to discuss this phenomenon in some generality. Consider 
a topological system supporting 27V Majorana zero-modes 
7i,...,2Ar. As in Eq. ( [98] ) we will (arbitrarily) combine the Ma- 
joranas into operators fj = (72^-1 + i72j)/2 whose occu- 
pation numbers rij = fjfj can be used to label the ground- 
state manifold. Suppose that one prepares this system into 
a ground state = |ni, . . . ,nAr) and then adiabatically 
swaps the positions of any two Majoranas. (We assume that 
the ground-state degeneracy is preserved throughout and that 
the initial and final Hamiltonians coincide.) The exchange 
statistics of the defects binding these zero-modes follows from 
the time evolution of |<l>^); for a nice discussion, see Refs. 12 
and|22l One source of this evolution is the dynamical phase 
e~i/o Eit)dt acquired by the wavefunction, where E{t) is 
the instantaneous ground- state energy during the interval T 
over which the interchange occurs. This factor is irrelevant 
for our purposes and will henceforth be ignored. More impor- 
tantly, the degeneracy together with the fractionalized nature 



of the zero-modes allows the system to end in a fundamen- 
tally different ground state \ ^f) from which it began. Proving 
this and the non-Abelian statistics that follows is nontrivial, 
and requires tracking the adiabatic evolution of the full many- 
body wavefunction — along with Berry matrices that can con- 
nect different ground states. Fortunately, one can deduce the 
final state | <!>/) (up to an overall phase) by addressing the dras- 
tically simpler problem of how the Majorana operators trans- 
form under the exchange. In all cases that we are aware of 
this procedure agrees with more rigorous approaches devel- 
oped, e.g., in Refs.[ini[Il|2ll[2ai32H3S| 

With this simplification in mind, we now review Ivanov'J^ 
remarkably accessible picture for non-Abelian statistics in a 
2D spinless p -\- ip superconductor (this discussion applies 
equally well to any of the experimental realizations in Sec. 
^V\. It suffices to analyze the configuration of Fig. fTSla) 



where four well-separated vortices trap Majorana zero-modes 
71,2,3,4. Suppose that we apply local pinning potentials to 
adiabatically braid the left two vortices clockwise as shown 
in the figure. This has two important consequences for the 
zero-modes — (i) 71 and 72 swap positions and (ii) 71 crosses 
a branch cut and acquires an additional minus sign. Thus the 
Majorana operators transform according to 71 —72 and 
72 ^ 7i • The unitary operator implementing this transforma- 
tion is U12 = (1 + 7i72)/V2 (that is, /7i27i/2^i2 = T72/i)- 
Similarly, clockwise exchange of neighboring vortices bind- 
ing 7j and 7^+1 sends — 7j+i and 7^+1 jj, which is 
generated by 



(l+7i7i+i)/v^- 



(111) 
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(For counterclockwise braids U~J_^-^ implements the transfor- 
mation.) Up to a common overall phase factor the ground- 
state wavefunctions |ni, 77,2) therefore evolve aJ^ |^i, ^2) 
Uj,j^i\ni,n2). 

For the three types of clockwise exchanges between neigh- 



boring vortices in Fig. 13 a), one explicitly finds 

|ni,n2) ^ /7i2|ni,n2) = e^^^^-^^^) Im, n2) (112) 
1^1,712) U23\ni,n2) (113) 

= -^[|ni,n2)+i(-l)"^|l-ni, 1-712)] 
|ni,n2)^/734|ni,n2) = 6^^(1-2^^) 



ni,n2). 



(114) 



Braiding 71^2 or 73^4 in a sense 'internally rotates' the ordi- 
nary fermion operators fj and produces nontrivial phase fac- 
tors in the states |n 1,712). More interestingly, braiding 72,3 
swaps 'half of /i with 'half of /2, resulting in a nontriv- 
ial rotation of |7ii, 77,2) within the ground-state manifold. To- 
gether these properties give rise to non-Abelian statistics of 
vortices: if one performs sequential exchanges, the final state 
depends on the order in which they are carried out. Math- 
ematically, this fascinating result follows from the nontrivial 
commutation relations satisfied by the operators in Eq. 
These conclusions generalize trivially to systems supporting 
arbitrarily many vortices. 

If one tries to extend this analysis to Majorana zero-modes 
arising in ID topological superconductors, two immediate 
problems arise. The first is that exchange statistics, non- 
Abelian or otherwise, is never well-defined in ID systems. 
As an example, suppose we attempt to adiabatically exchange 
73 and 74 in the setup of Fig. |6jd) by moving 73 rightward 
and 74 leftward. Clearly these zero-modes eventually overlap 
and split the ground-state degeneracy; furthermore, whether 
or not one actually performed an exchange when the system 
returns to its original configuration is completely ambiguous. 
Moving away from strict one-dimensionality by fabricating 
networks of ID wireJSlor 2D topological insulator edges^^^ 
circumvents this problem in a conceptually straightforward 
manner. Figures [T3jb) and (c) illustrate examples that allow 
Majorana modes to be meaningfully exchanged. (As usual 
green denotes topological regions while dashed lines are triv- 
ial. Note also that the structure of such networks is rather 
arbitrary — they can even form three-dimensional lattices. 
The zero-modes can be adiabatically transported in these se- 



tups by applying gate voltage^^^*^ or supercurrents^^^ to shift 
the domain wall locations as desired. In this way one can 



turn 
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exchange 71 and 72 in Fig. 13 ^b) by executing a 'three-point- 



first moving 71 to the center of the vertical line, 
then moving 72 all the way leftward, and finally moving 71 
up and to the right. Similar ideas allow one to exchange 
the Majoranas 73 and 74 belonging to different topological 
segments, and can also be adapted to the corner junction of 
Fig. 13 c). Interestingly, novel methods of effectively imple- 



menting exchanges without physically transporting Majoranas 
(as in measurement-only topolog ical qua ntum computation"^^) 
have also been proposed recently.ESESl 

Although the exchange of Majorana zero-modes becomes 
well-defined in these networks, a second, much subtler prob- 



lem appears in this context. Namely, the vortices that are cru- 
cial for establishing non-Abelian statistics in 2D p + ip su- 
perconductors are entirely absent here — so does non-Abelian 
statistics still emerge? Fortunately several studies have shown 
that it does, and that the Majorana zero-modes in fact trans- 
form under exchange exactly as in the p -\- ip case despite 
the lack of vortices. ^^ ^^ To provide a rough flavor for how 
this arises, consider the physical situation where the network 
of Fig. \T3[h) arises from spin-orbit-coupled wires adjacent to 
an 5-wave superconductor with uniform phase (j). Suppose 
that we exchange 71^2 as described above. Although the 5- 
wave superconductor exhibits a uniform phase, as the Majo- 
ranas traverse the network the phase of the effective p-wave 
pair field that they experience does in fact vary. This variation 
causes one of the Majoranas to acquire a minus sign arising 
from a branch cut, precisely as in Ivanov's construction.^^ For 
more details and complementary perspectives on this interest- 
ing problem see Refs. [35H38l 

By virtue of non-Abelian statistics, braiding Majorana zero- 
modes in 2D p ^ ip superconductors and ID networks al- 
lows one to perform topological quantum information pro- 
cessing that is in principle immune from decoherence.^^ Equa- 



tions \\\2) through illustrate a concrete example of 

the protected qubit processing effected in this manner. Un- 
fortunately, such qubit rotations are too restrictive to per- 
mit universal quantum computation; two additional processes 
are needed. The first is a 7r/8 phase gate that introduces 
phase factors e^*^/^ depending on the occupation number 
corresponding to a given pair of Majoranas. The second is 
the ability to read out the eigenvalue of the product of four 
Majoranas, jijjjklh without measuring that of individual 
pairs. While these processes can introduce errors this does 
not mean that the topologically protected braiding operations 
are without merit. In fact the unprotected part of the compu- 
tation enjoys a dramatically higher error thres hold compared 
to conventional quantum computing schemes. '^^'^ Many in- 
genious proposals already exist for supplement braiding with 
the operations required to perform universal quantum compu- 
tation, both in arch itectures based on ID and 2D topological 
phases 1 8 | 255 | 257j [267 blueprints for a Majorana-based 
quantum computer are therefore already in place; we simply 
need to begin assembling the hardware. 



VII. OUTLOOK 

The possibility of observing Majorana fermions in con- 
densed matter systems now appears tantalizingly close. In 



Sees. [Ill] and |IV] we saw that the number of realistic pro- 
posals that now exist is rather immense, and most involve 
heterostructures with garden- variety 5- wave superconductors. 
One point worth emphasizing here is that the theory for 
these proximity-induced topological superconductors centers 
largely around non-interacting electron models — apart from, 
of course, attractive interactions implicitly invoked in the par- 
ent superconductors. We view this simplicity as an enormous 
virtue that affords theorists a degree of predictive power for 
experiments that is not often encountered in the quest for ex- 
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otic phases of matter. (One can argue that the rapid progress 
in the field of topological insulators arose for similar reasons.) 
Given this fact, along with the accessibility of the building 
blocks comprising the heterostructures, it is no surprise that 
several of these new proposals have already inspired broad 
experimental efforts in leading laboratories worldwide. If 
these efforts continue unabated we believe the question is not 
whether Majorana modes will be definitively identified, but 
rather when, in what system, and with what measurement. 

Amongst the numerous proposals reviewed the original 
topo logic al-insulator-based devices introduced by Fu and 
Kane3^E^(see also Ref. ll62t remain in some ways ideal. The 
elegance of the theories is hard to match, and from a prac- 
tical standpoint the prospect of obtaining topological phases 
protected by gaps that are immune to disordei^"^ and lim- 
ited only by that of the parent superconductoi^^ is exceed- 
ingly attractive. We believe strongly that experiments in this 
direction should be ardently pursued using both 2D and 3D 
topological insulators even if another system 'wins the race' 
for Majorana; future applications may rest on such excep- 
tional properties. Proposals employing conventional spin- 
orbit-coupled ID wires, first pioneered in Refs. lll5l andlll61 
also stand out given the comparative maturity of semicon- 
ductor technology as well as the simplicity and tunability 
of the required architectures. These systems are well-poised 
to experimentally realize the predicted magnetic-field-driven 
topological phase in the near future. ^^^ ^^^ More broadly, we 
hope that experimentalists will push many more of the newly 
proposed directions while theorists continue to conjure up 
new and improved Majorana platforms. There may just be 
a sleeper in the mix, or perhaps the ideal direction is even yet 
to be introduced. 

The first unambiguous sighting of Majorana fermions in 
condensed matter would provide a landmark event in physics. 
We hope to have made a compelling case in Sec. |V] that nu- 
merous smoking-gun detection methods are now available to 
make this identification definitive. It should be stressed that 
this initial observation will only herald the beginning of what 
is likely to be a long, fruitful subfield. The realization of 
exotic physics such as fractional Josephson effects and non- 
Abelian statistics, as well as applications from topological 
quantum memory to universal quantum computation are truly 
fascinating goals that will keep physicists occupied for many 
years to come. And as always there are bound to be many 
surprises along the way. 

To conclude we will briefly highlight some interesting fu- 
ture directions and open questions. One intriguing alterna- 
tive route to engineering topological phases involves period- 
ically driving a system that would otherwise be trivial. Us- 
ing this mechani sm prop osals for artificially generating topo- 
logical insulator JSmir ^j^^ s^perfluids supporting Ma- 
jorana fermion^^ have recently been put forth. The pos- 
sibility of moving away from static systems as a means of 
generating non-Abelian phases opens new avenues that war- 
rant further exploration. There are also many other promis- 
ing routes to Majorana fermions that we have not touched on 
here. Among the most interesting is the potential realization 
of Kitaev's 2D honeycomb model^^^ in a certain class of mag- 



netic insulators^jt^. Looking forward, it is worth exploring 
whether the connection identified by Read and Green be- 
tween the Moore-Read state^^ and a spinless p-\-ip supercon- 
ductor can be adapted to still more exotic fractional quantum 
Hall states. As noted in the introduction it is this remark- 
able correspondence that led to the realization that a weakly 
correlated superconductor can harbor non-Abelian statistics. 
Establishing a similar correspondence between quantum Hall 
phases supporting even richer non-Abelian anyons and phases 
exhibited by less strongly interacting 2D systems could open 
entirely new directions in the pursuit of topological quantum 
computation. If this can be achieved, might there exist related 
ID systems supporting these richer non-Abelian anyons, in 
the same way that Majorana modes can appear in either ID or 
2D topological superconductors? 
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Appendix A: Derivation of the effective action for ID and 2D 
systems coupled to a bulk s-wave superconductor 

Consider a d-dimensional system of electrons (with d = 
1 or 2) proximate to a bulk 5-wave superconductor. Let the 
Hamiltonian for this structure he H = Hd -\- Hsc + ^r, 
where 



-I 



(Al) 



describes the (i-dimensional system, 

J7^[^kriir]i^ + A5c(/?tk^i-k + i^.c.)], (A2) 

models the 5-wave superconductor, with ek = k'^/{2msc) — 
jisc the superconductor's kinetic energy, and 

^(^/'k,7?k + i^.c.) (A3) 

the term which incorporates electron tunneling between the 



two subsystems. In Eq. (A3), k^^ = kx if d = 1 while k^^ 
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{kx^ky) if d = 2. Our goal is to obtain an effective action 
for the d-dimensional system with the gapped superconductor 
degrees of freedom integrated out. To achieve this it is useful 
to first perform a unitary transformation which diagonalizes 
Hsc- 



,(A4) 



where = y^e^ + are the quasiparticle energies for the 
superconductor. In this new basis one obtains 



Hsc = J j^Ek[x[i^Xi\c + X2kX2k] 

^^;^[XikKV^jk^ + ^feV^i-kJ 



(A5) 



(A6) 



It is now straightforward to write down the Euclidean path 
integral corresponding to H and then integrate out the quasi- 
particle operators xi,2- This yields an effective action S^f^ = 
Sd + SS, with 



5S 



/duj (i^k 
^(^%,-)'^k%,o.) (A7) 

duJ d^k r A A /T N r / / rr 1 

|A5cA(k,Cj)[V;^(k,a;)V;;(_k -u;) + H.C.\ 



27T (27r)^ 

+ [-iujX{k,uj) - (5/i(k,a;)]V^[j^^^^%,u;)}. (A8) 



One sees here that the superconductor renormalizes the chem- 
ical potential for the (i-dimensional system (in a weakly 
frequency- and momentum-dependent fashion) through 



^/i(k, cj); this correction is unimportant, however, and will be 
henceforth neglected. The essential physics associated with 
the hybridization is encoded in the function A(k,cj) appear- 
ing in SS. Let us focus for the moment on d = 1 where this is 
given by 



(A9) 

Typically we will be concerned with one-dimensional systems 
of rather low density so that kl/{2msc) ^ Msc for the impor- 
tant values of kx ; in this case the dependence of A on kx can 
be safely igno red. Maki ng the further reasonable assumption 
that jiisc ^ aA^^^-^AI^ over the relevant frequencies, A eval- 
uates to the following simple expression: 



X{oj) 



A2 



(AlO) 



where p2D = i^sc/ (27r) is the density of states (per spin) for a 
two-dimensional system with effective mass rrisc- For a d = 2 
dimensional system \ {kx-,ky^uj) follows from the right side of 
Eq. \A9) integrated only over kz. Under similar assumptions 
made for d = 1 A is again approximately momentum inde- 



pendent and given by Eq. (AlO) with p2D replaced by the ID 
density of states piD = \/'^scl i^^sc]- Upon defining a 
quasiparticle weight Z{uj) = [1 + A(a;)]~^, the effective ac- 
tion can be expressed in the desired form quoted in the main 
text: 



/duj d^k. 
^^^:^^"^M{V^(k,u;)[-^^ + ^M'^k]%,u;) 

+ A,e[l - Z(cj)][V^^(k,a;)V^;(_k,-a;) +i^.C.]}. (All) 
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